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Abstract 



In this paper, we show that the adaptive multidimensional increment ratio estimator of the long range 
memory parameter defined in Bardet and Dola (2012) satisfies a central limit theorem (CLT in the sequel) 
for a large semiparametric class of Gaussian fractionally integrated processes with memory parameter 
d € (—0.5, 1.25). Since the asymptotic variance of this CLT can be computed, tests of stationarity or 
nonstationarity distinguishing the assumptions d < 0.5 and d > 0.5 are constructed. These tests are 
also consistent tests of unit root. Simulations done on a large benchmark of short memory, long memory 
and non stationary processes show the accuracy of the tests with respect to other usual stationarity or 
nonstationarity tests (KPSS, LMC, ADF and PP tests). Finally, the estimator and tests are applied to 
. log-returns of famous economic data and to their absolute value power laws. 

. Keywords: Gaussian fractionally integrated processes; Adaptive semiparametric estimators of the meme- 

ory parameter; test of long-memory; stationarity test; unit root test. 
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1 Introduction 

Consider the set I{d) of fractionally integrated time series X = {Xk)kez for —0.5 < d < 1.5 by: 

Assumption I{d): X = {Xt)tez is a time series if there exists a continuous function f* : [— 7r,7r] — > [0, oo[ 
satisfying: 

1. if —0.5 < d < 0.5, X is a stationary process having a spectral density f satisfying 
/(A) = |A|-2'*/*(A) for all A e (-7r,0)U (0,7r), with /*(0) > 0. (1.1) 

2. i/0.5 < d < 1.5, U = {Ut)tez ~ Xt — Xt-i is a stationary process having a spectral density f satisfying 
/(A) = |A|2"^''/*(A) for ah Ae (-7r,0)U (0,7r), with /*(0) > 0. (1.2) 

The case d E (0, 0.5) is the case of long-memory processes, while short-memory processes are considered 
when —0.5 < d < and nonstationary processes when d > 0.5. ARFlMA{p, d, q) processes (which are linear 
processes) or fractional Gaussian noises (with parameter H = d + 1/2 E (0, 1)) are famous examples of pro- 
cesses satisfying Assumption I{d). The purpose of this paper is twofold: firstly, we establish the consistency 
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of an adaptive semiparametric estimator of d for any d e (—0.5, 1.25). Secondly, we use this estimator for 
building new semiparametric stationary tests. 

Numerous articles have been devoted to estimate d in the case d € (—0.5,0.5). The books of Beran (1994) 
or Doukhan et al. (2003) provide large surveys of such parametric (mainly maximum likelihood or Whittle 
estimators) or semiparametric estimators (mainly local Whittle, log-periodogram or wavelet based estima- 
tors). Here we will restrict our discussion to the case of semiparametric estimators that are best suited to 
address the general case of processes satisfying Assumption I{d). Even if first versions of local Whittle, log- 
periodogramm and wavelet based estimators (see for instance Robinson, 1995a and 1995b, Abry and Veitch, 
1998) are only considered in the case d < 0.5, new extensions have been provided for also estimating d when 
d > 0.5 (see for instance Hurvich and Ray, 1995, Velasco, 1999a, Velasco and Robinson, 2000, Moulines and 
Soulier, 2003, Phillips and Shimotsu, 2005, Giraitis et ai, 2003 or Moulines et ai, 2007,). Moreover, adap- 
tive versions of these estimators have also been defined for avoiding any trimming or bandwidth parameters 
generally required by these methods (see for instance Moulines and Soulier, 2003, or Veitch et at, 2003, 
or Bardet et ai, 2008). However there still no exists an adaptive estimator of d satisfying a central limit 
theorem (for providing confidence intervals or tests) and valid for d < 0.5 but also for d > 0.5. This is the 
first objective of this paper and it will be achieved using multidimensional Increment Ratio (IR) statistics. 
Indeed, Surgailis et al. (2008) first defined the statistic IRn (see its definition in (|2.3p ') from an observed 
trajectory {Xi, . . . , X^). Its asymptotic behavior is studied and a central limit theorem (CLT in the sequel) 
is established for d e (—0.5,0.5) U (0.5,1.25) inducing a CLT. Therefore, the estimator dN — Aq^{IRn), 
where d i— >■ Aq (d) is a smooth and increasing function, is a consistent estimator of d satisfying also a CLT (see 
more details below). However this new estimator was not totally satisfying: firstly, it requires the knowledge 
of the second order behavior of the spectral density that is clearly unknown in practice. Secondly, its nu- 
merical accuracy is interesting but clearly less than the one of local Whittle or log-periodogram estimators. 
As a consequence, in Bardet and Dola (2012), we built an adaptive multidimensional IR estimator d)^ 
(see its definition in p.2p ) answering to both these points but only for —0.5 < d < 0.5. This is an adaptive 
semiparametric estimator of d and its numerical performances are often better than the ones of local Whittle 
or log-periodogram estimators. 

Here we extend this preliminary work to the case 0.5 < d < 1.25. Hence we obtain a CLT satisfied by d)^ 
for all d G (—0.5, 1.25) with an explicit asymptotic variance depending only on d and this notably allows to 
obtain confidence intervals. The case d = 0.5 is now studied and this offers new interesting perspectives: our 
adaptive estimator can be used for building a stationarity (or nonstationarity) test since 0.5 is the "border 
number" between stationarity and nonstationarity. 

There exist several famous stationarity (or nonstationarity) tests. For stationarity tests we may cite the 
KPSS (Kwiotowski, Phillips, Schmidt, Shin) test (see for instance Hamilton, 1994, p. 514) and LMC test 
(see Leybourne and McCabe, 2000). For nonstationarity tests we may cite the Augmented Dickey- Fuller 
test (ADF test in the sequel, see Hamilton, 1994, p. 516-528) and the Philipps and Perron test (PP test in 
the sequel, see for instance Elder, 2001, p. 137-146). All these tests are unit root tests, i.e. and roughly 
speaking, semiparametric tests based on the model Xt = pXt-i + £t with \p\ < 1. A test about d = 0.5 for a 
process satisfying Assumption I{d) is therefore a refinement of a basic unit root test since the case p = 1 is a 
particular case of /(I) and the case |/5| < 1 a particular case of /(O). Thus, a stationarity (or nonstationarity 
test) based on the estimator of d provides a more sensible test than usual unit root tests. 
This principle of stationarity test linked to d was also already investigated in many articles. We can no- 
tably cite Robinson (1994), Tanaka (1999), Ling and Li (2001), Ling (2003) or Nielsen (2004). However, aU 
these papers provide parametric tests, with a specified model (for instance ARFIMA or ARFIMA-GARCH 
processes). To our knowledge, such tests were not yet proposed and studied in a general semiparametric 
framework. 

Here we slightly restrict the general class I{d) to the Gaussian semiparametric class IG{d, /3) defined below 
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(see the beginning of Section For processes belonging to this class, we construct a new stationarity 
test Sjy which accepts the stationarity assumption when dj^ < 0.5 + s with s a threshold depending on 
the type I error test and N, while the new nonstationarity test Tn accepts the nonstationarity assumption 
when rf^^'' > 0.5 — s. Note that < s' also provides a test for deciding between short and long range 

dependency, as this is done by the V/S test (see details in Giraitis et ai, 2003) 

In Section [5l numerous simulations are realized on several models of time series (short and long mem- 
ory processes). 

First, the new multidimensional IR estimator d^^'' is compared to the most efficient and famous semipara- 
metric estimators for d G [—0.4, 1.2]; the performances of djy are extremely convincing and globally better 
than the other estimators. 

Secondly, the new stationarity Sn and nonstationarity Tn tests are compared on the same benchmark of 
processes to the most famous unit root tests (KPSS and LMC, ADF and PP tests). And the results are quite 
surprising: even on AR[1] or ARIMA[1, 1, 0] processes, multidimensional IR Sn and T/v tests provide better 
results than KPSS and LMC stationarity tests which are nevertheless constructed for such processes. Note 
however that ADF and PP tests provide results slightly better than Sn and Tjv tests for these processes. 
For long- memory processes (such as ARFIMA processes), the results are clear: Sn and Tn tests are efficient 
tests of (non)stationarity while KPSS, LMC, ADF and PP tests are not relevant at all. 
Finally, we studied the stationarity and long range dependency properties of Econometric data. We chose to 
apply estimators and tests to the log-returns of daily closing value of 5 classical Stocks and Exchange Rate 
Markets. After cutting the series in 3 stages using an algorithm of change detection, we found again this well 
known result: the log- returns are stationary and short memory processes while absolute values or powers 
of absolute values of log-returns are generally stationary and long memory processes. Classical stationarity 
or nonstationarity tests are not able to lead to such conclusions. We also remarked that these time series 
during the "last" (and third) stages (after 1997 for almost all) are generally closer to nonstationary processes 
than during the previous stages with a long memory parameter close to 0.5. 

The forthcoming Section [2] is devoted to the definition and asymptotic behavior of the adaptive multidi- 
mensional IR estimator of d. The stationarity and nonstationarity tests are presented in Section |4] while 
Section [S] provides the results of simulations and application on econometric data. Finally Section [S] contains 
the proofs of main results. 

2 The multidimensional increment ratio statistic 

In this paper we consider a semiparametric class IG{d, (3): for < d < 1.5 and /3 > define: 

Assumption IG{d,l3): X = {Xt)t£Z is a Gaussian time series such that there exist e > 0, cq > 0, 
Cq > and ci G R satisfying: 

1. if d < 0.5, X is a stationary process having a spectral density f satisfying for all A G (— 7r,0) U (0,7r) 

/(A) = co|Ar2'^ + ci|Ar2'i+/3 + o(|A|-2W.) and |/'(A)| < cj, A-^-i-i. (2.1) 

2. i/0.5 < d < 1.5, U = {Ut)te'i, = Xt-i is a stationary process having a spectral density f satisfying 
for all A e (-7r,0)U (0, tt) 

/(A)==co|A|2-2'i + ci|A|2-2'i+'3 + 0(|A|2-2W.) ^ud | /'(A) | < c[, A-^'i+i . (2.2) 

Note that Assumption IG{d, /?) is a particular (but still general) case of the more usual set I{d) of fractionally 
integrated processes defined above. 
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Remark 1. We considered here only Gaussian processes. In Surgailis et al. (2008) and Bardet and Dola 
(2012), simulations exhibited that the obtained limit theorems should be also valid for linear processes. How- 
ever a theoretical proof of such result would require limit theorems for functionals of multidimensional linear 
processes difficult to be established. 



In this section, under Assumption IG{d,j3), we establish central limit theorems which extend to the case 
d e [0.5, 1.25) those already obtained in Bardet and Dola (2012) for d € (-0.5,0.5). 

Let X = {Xk)kQn be a process satisfying Assumption IG{d, P) and [Xi,- ■ ■ , Xn) be a path of X. For any 
ieN* define 

k+i k+e k+2e k+2i 

rp (e\ — 1 \ ^ t=k+i t=fc+i t=k+i+i t=k+t+i „x 

^^>-~N-U k+l k+t k+2i k+2e • ''^■'^) 



fe=0 / \ - ^ \ ^ 



Xt+,- ^ Xt) +\i >^ Xt+e- >^ Xt) 



t=fe+l t=fe+l t=k+l+\ t=k+t+\ 

The statistic IRn was first defined in Surgailis et al. (2008) as a way to estimate the memory parameter. 
In Bardet and Surgailis (2011) a simple version of IR-statistic was also introduced to measure the roughness 
of continuous time processes. The main interest of such a statistic is to be very robust to additional or 
multiplicative trends. 

As in Bardet and Dola (2012), let mj = j m, j = 1, - ■ ■ ,p with p e N* and m e N*, and define the random 

vector {IRN(mj))i<j<p- In the sequel we naturally extend the results obtained for m €E N* to m S (0, oo) 
by the convention: {IR]\r{j mj)i<j<p = (/i?jv(i [m]))i<j<p (which change nothing to the asymptotic results). 

For H £ (0,1), let Bh = {BH{t))teR be a standard fractional Brownian motion, i.e. a centered Gaus- 
sian process having stationary increments and such as Cov(Sji-(f) , Bh{s)) = | + ~ I* ~ •^1^^)- 
Now for d e (-0.5,0.5) U (0.5, 1.25) and j e N*, define for r e M, 



^ ^\ [\Bd-Ur + s+j)-B,_Ur + s))ds if d S (0.5, 1.25) 

1^ ^^^^^^ ^ _ (Bd+o.5(r + 2j) - 2 Bd+o.5{T + j) + Bd+o.5{T)) if d e (-0.5, 05) 

Note that (Z^"'^(t))^ is a stationary Gaussian process for any d € (—0.5,0.5) U (0.5,1.25). When d = 
0.5, such definition does not hold. However for p e N*, a definition of the multidimensional process 
{Zq^I{t), . . . , Zq^I{t)) can be obtained by continuous extension when d 0.5 of the covariance of 
{Z^'\0),Z^J\t)). Therefore for t e K, we define {Z^^I{t), . . . , Z^^},{t)) ^ as a centered stationary Gaus- 
sian multidimensional process satisfying: 

Cov(Z« (0), 4^)(r)) := (^-h{T + i-j) + h{T + i) + h{T - j) - /i(t)) for r e M, 

with h{x) = |(|a; — l|^log|a; — 1| -|- |a; -|- lplog|a; -l- 1| — 2|a;|^ log |a;|) for a; e M, using the convention 
X logO = 0. Now, we establish a multidimensional central limit theorem satisfied by (/-Rjv(j 'Ti))i<j<p for 
all d e (-0.5, 1.25): 

Proposition 1. Assume that Assumption IG{d,P) holds with —0.5 < d < 1.25 and (3 > 0. Then 

M(lRNijm)-E[lRMUm)]) A ^(0,rp(d)) (2.5) 

with Tp{d) = {o-ij{d))i<ij<p where for i,j e {1, . . . ,p}, 
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The proof of this proposition as well as all the other proofs is given in Section [51 As numerical experiments 
seem to show, we will assume in the sequel that Tp{d) is a definite positive matrix for all d G (—0.5, 1.25). 
Now, this central limit theorem can be used for estimating d. To begin with. 

Property 2.1. Let X satisfying Assumption IG{d,/3) with 0.5 < d < 1.5 and < /3 < 2. Then, there exists 
a non-vanishing constant K{d,l3) depending only on d and /3 such that for m large enough, 



E[lRNim)] 



Ao{d) + K{d, P) X (1 + o(l)) ifl3<l + 2d 

Ao(d) +i^(0.5,/3) xm-^\ogm{l + o{l)) if ^ ^ 2 and d = 0.5 



Ad+1.5 _ n<i+0.5 _ j 

for 0.5 < d < 1.5 



wtth Ao(d) := A{p{d)) where p{d) :^ { g (2-7) 

TT, — 7777 — 2 for d = 0.5 

81og(2) 



A(r):=-arctanJi±^ + i,/i±^log(-^) for \r\ < 1. (2.8) 

TT y 1 — r TT y 1 — r 1 + r 



Therefore by choosing m and A'' such as (^y/N/mjrn^^ log rn — > when m,N ^ oo, the term E[/i?(jm)] 
can be replaced by Ao(d) in Proposition [TJ Then, using the Delta-method with the function {xi)i<i<p 
{AQ^{xi))i<i<p (the function d € (—0.5, 1.5) — > Ao((i) is a C°° increasing function), we obtain: 

Theorem 1. Let dj^Ufn) := A^^ (^IR^ij rn)^ for I < j < p. Assume that Assumption IG{d,(3) holds with 
0.5 < d < 1.25 and < ^ < 2. Then if m ^ C m with C > and (1 + 2/3)-i < a < 1, 



\[^(Mjm)-d) A M(o,{A',id))-'rp{d) 

V mV / i<j<p N^oc \ 



(2.9) 



This result is an extension to the case 0.5 < d < 1.25 from the case —0.5 < d < 0.5 already obtained in 
Bardet and Dola (2012). Note that the consistency of dAr(jm) is ensured when 1.25 < d < 1.5 but the 
previous CLT does not hold (the asymptotic variance of y^dAr(jm) diverges to oo when d — > 1.25, see 
Surgailis et al, 2008). 
Now define 

EAr(m) :- {A'a{dN{m))-^rp{dN{m)). (2.10) 
The function d G (—0.5, 1.5) i-^ cr(d)/A'(d) is C°° and therefore, under assumptions of TheoremfH 



S^(m) A {A'oid))-^Tp{d). 



Thus, a pseudo-generalized least square estimation (LSE) of d ican be defined by 

dAr(TO) := (j^(EAr(m)) Vp) ^ (S7v(m)) ^ (rfjv(n^i)) i<i<p 

with Jp := (l)i<j<p and denoting its transpose. From Gauss-Markov Theorem, the asymptotic variance 
of dN{m) is smaller than the one of dj\[{jfn), j — 1, . . . ,p. Hence, we obtain under the assumptions of 
Theorem [TJ 



(d^(TO)-d) aa(o, A;,(d)-2(j;rpi(d)jp) (2.11) 
3 The adaptive version of the estimator 

Theorem [1] and CLT (|2.1ip require the knowledge of (3 to be applied. But in practice /? is unknown. The 
procedure defined in Bardet et al. (2008) or Bardet and Dola (2012) can be used for obtaining a data-driven 



5 



selection of an optimal sequence (toat) derived from an estimation of 13. Since the case d G (—0.5,0.5) was 
studied in Bardet and Dola (2012) we consider here d G [0.5, 1.25) and for a G (0, 1), define 

QNia,d) := {Mj N'^) -dNiN")y^^^^^{^MiN'^)y\dMij N'^) -MN'^))^^^^^^ (3.1) 

which corresponds to the sum of the pseudo- generalized squared distance between the points {dM{j and 
PGLS estimate of d. Note that by the previous convention, iInU iV") = dN{j [N"]) and dAr(iV") = djv([^"])- 
Then QN^ct) can be minimized on a discretization of (0, 1) and define: 



aN 



Argmin„g_4^QAr(Q;) with An ^ ^ 



_^ 3_ \ogN/p\_^ 

\osN ' logA^ losN /■ 



Remark 2. The choice of the set of discretization An is implied by our proof of convergence of q.n to 
a* . If the interval (0, 1) is stepped in N'^ points, with c > 0, the used proof cannot attest this convergence. 
However log N may he replaced in the previous expression of An by any negligible function of N compared 
to functions N'^ with c > (for instance, (logiV)" or alogN with a > can be used). 

From the central limit theorem (12.91) one deduces the following limit theorem: 

Proposition 2. Assume that Assumption IG[d,(3) holds with 0.5 < d < 1.25 and < /3 < 2. Then, 

^ V . 1 

aN — > 01 = 



N^oo (1 + 2/J)' 

Finally define 

3„ -.v, ~ , log log iV 

TOAr:=iV" with aN:^aN + -, ttt- ■ 7^- 

(p- 2)(1 - ajv) logiV 

and the estimator 

^^""^ := dN{mN) = djv(iV""). (3.2) 

(the definition and use of ctN instead of Sat are explained just before Theorem 2 in Bardet and Dola, 2012). 
The following theorem provides the asymptotic behavior of the estimator dj^ ■ 

Theorem 2. Under assumptions of Proposition\^ 

-d) AA(o; A^(d)-2 (j;r^-i(d)J,)-^). (3.3) 



M w . 2(1 + 3/3) NTTW 

Moreover, Vp > — — , — ■ \d'^ ' - d\ — > 0. 

(p—2)p (\0g NjP ' ' N^oo 



The convergence rate of d^^"* is the same (up to a multiplicative logarithm factor) than the one of minimax 
estimator of d in this semiparametric frame (see Giraitis et al., 1997). The supplementary advantage of 
d^]^^^ with respected to other adaptive estimators of d (see for instance Moulines and Soulier, 2003, for 

— i T 

an overview about frequency domain estimators of d) is the central limit theorem (13.31) satisfied by dj^ . 
Moreover d)^ can be used for d e (—0.5,1.25), i.e. as well for stationary and non-stationary processes, 
without modifications in its definition. Both this advantages allow to defined a stationarity test based on 



4 New stationarity and nonstationarity tests 

Assume that {Xi,...,Xn) is an observed trajectory of a process X ~ {Xk)kez- We define here new 
stationarity and nonstationarity tests for X based on (J]^^^ . 
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4.1 A new stationarity test 

There exist many stationarity and nonstationarity test. The most famous stationarity tests are certainly the 
following miit root tests; 

• The KPSS (Kwiotowski, Phillips, Schmidt, Shin) test (see for instance Hamilton, 1994, p. 514); 

• The LMC (Leybourne, McCabe) test which is a generalization of the KPSS test (sec for instance 
Leybourne and McCabe, 1994 and 1999) 

Here, we consider the following problem of test: 

• Hq: X is a stationary process. 

• Hi: X is not a stationary process. 

We use a test based on (i^^"* for deciding between these hypothesis. Hence from the previous CLT 13.31 and 
with a type I error a, define 

Sn ■■= l4^^«)>o.5+o-p(o.5)gi_„JV(°«-i)/=^' (4-1) 

where crp(0.5) = (A;,(0.5)-2 ( 1(0.5) Jp) ] (see ^) and qi^a is the (1 — a) quantile of a standard 

Gaussian random variable J\f{0, 1). 

Then we define the following rules of decision: 

• Hq is accepted (or equivalently X is a stationary process) when Sn = 0. 

• Hq is rejected (or equivalently X is a nonstationary process) when Sn — 1- 
From previous results, it is clear that: 

Property 1. If X satisfies Assumption IG{d,/3) with d < 0.5 and /3 e (0,2] (stationary process), the 
asymptotic type I error of the test Sn is a and the test power tends to 1 if X Cz IG{d, j3) with d > 0.5 and 
/?€ (0,2]. 

Moreover, this test can be used as a unit root test. Indeed, define the following typical problem of unit 
root test. Let Xt — at + b + et, with (a, 6) G M^, and et an ARIMA(p, d, q) with d = Q oi d = 1. Then, a 
(simplified) problem of a unit root test is to decide between: 

• Hq^: d — and (e^) is a stationary ARMA(p, q) process. 

• Hi^: d = 1 and {et — £t-i)t is a stationary ARMA(p, q) process. 
Then, 

Property 2. Under assumption Hq^, the type I error of this unit root test problem using Sn decreases to 
when Af — >■ oo and the test power tends to 1. 

4.2 A new nonstationarity test 

Famous unit root tests are more often nonstationarity test. For instance, between the most famous tests, 

• The Augmented Dickey-Fuller test (see Hamilton, 1994, p. 516-528 for details); 

• The Philipps and Perron test (a generalization of the ADF test with more lags) (see for instance Elder, 
2001, p. 137-146). 

~^ IR) — 

Using the statistic djy we propose a new nonstationarity test Tn for deciding between: 
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• Hq. X is a nonstationary process when — 1; 

• H[: X is a stationary process when = 0, 
where 

Tjv := ldU«)<o.5-,.p(0.5)gi-„Ar(°iv-i)/2- (4-2) 

Then as previously 

Property 3. If X satisfies Assumption IG{d,j3) with d > 0.5 and (3 G (0,2] (nonstationary process), the 
asymptotic type I error of the test Tn is a and the test power tends to 1 if X Cz IG{d,l3) with d < 0.5 and 

/Je (0,2]. 

As previously, this test can also be used as a unit root test where Xt = at + b + St, with (a, b) e K^, and 
Et an ARIMA(p, d, q) with d = or = 1. We consider here a "second" simplified problem of unit root test 
which is to decide between: 

• Hq^ : d — 1 and {et — £t-i)t is a stationary ARMA(p, q) process. 

• : d = and {et)t is a stationary ARMA(p, q) process.. 
Then, 

Property 4. Under assumption Hq^ , the type I error of the unit root test problem using decreases to 
when N ^ oo and the test power tends to 1. 

5 Results of simulations and application to Econometric and Fi- 
nancial data 

5.1 Numerical procedure for computing the estimator and tests 

First of all, softwares used in this Section are available on |http: //samm.univ-parisl . f r/-Jean-Marc-Bardet 

with a free access on (in Matlab language) as well as classical estimators or tests. 

The concrete procedure for applying our MIR-test of stationarity is the following: 

1. using additional simulations (realized on ARM A, ARFIMA, FGN processes and not presented here 
for avoiding too much expansions) , we have observed that the value of the parameter p is not really 
important with respect to the accuracy of the test (less than 10% on the value of d^^''). However, for 
optimizing our procedure we chose p as a stepwise function of n: 

p — 5 X 1{„<120} + 10 X l{i20<n<800} + 15 X l{800<n<10000} + 20 X 1{„>10000} 

and 0-5(0.5) ~ 0.9082; f7io(0.5) ~ 0.8289; 015(0.5) ~ 0.8016 and o-2o(0.5) ~ 0.7861. 

2. then using the computation of m^r presented in Section [3l the adaptive estimator 3]^^^ (defined in 
(j3.2p ) and the test statistics Sn (defined in (|4.ip ) and (defined in ()4.2p ) are computed. 

5.2 Monte-Carlo experiments on several time series 

In the sequel the results are obtained from 300 generated independent samples of each process defined below. 
The concrete procedures of generation of these processes are obtained from the circulant matrix method, as 
detailed in Doukhan et at (2003). The simulations are realized for different values of d and TV and processes 
which satisfy Assumption IG{d,j3): 
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1. the usual ARIMA(p, d, q) processes with respectively c? = or o? = 1 and an innovation process which 
is a Gaussian white noise. Such processes satisfy Assumption IG{0, 2) or IG{1, 2) holds (respectively); 

2. the ARFIMA(p, d, q) processes with parameter d such that d e (—0.5, 1.25) and an innovation process 
which is a Gaussian white noise. Such ARFIMA(p, d, q) processes satisfy Assumption IG{d, 2) (note 
that ARIMA processes are particular cases of ARFIMA processes). 

3. the Gaussian stationary processes X'^'^'^\ such as its spectral density is 

/3W = ^(1 + ^'') for Ae [-^,0)U(0,^], (5.1) 

with d € (—0.5,1.5) and f3 G (0,oo). Therefore the spectral density /a implies that Assumption 
IG{d, 13) holds. 

Comparison of d!"^^'^ with other semiparametric estimators of d 

^ I R) 

Here we first compare the performance of the adaptive MIR-estimator dpf with other famous semipara- 
metric estimators of d: 

• dH'j^^^^ is the adaptive global log-periodogram estimator introduced by Moulines and Soulier (2003), 
also called FEXP estimator, with bias- variance balance parameter k = 2. Such an estimator was shown 
to be consistent for d e] — 0.5, 1.25]. 

• is the local Whittle estimator introduced by Robinson (1995b). The trimming parameter is 
chosen as m = A^/30 (this is not an adaptive estimator). This estimator was shown to be consistent 
for d e] ~ 0.5,1] and not consistent for d > 1 (see Phillips and Shimotsu, 2004), with normality for 
d e] — 0.5,0.75[ (see Velasco, 1999b), consistency and normality for d e] — 0.5, 1.5[ with tapering only 
(see Hurvich and Chen, 2000). 

• (i^'^^^' is an adaptive wavelet based estimator introduced in Bardet et al. (2008) using a Lemarie- 
Meyer type wavelet (another similar choice could be the adaptive wavelet estimator introduced in Veitch 
et al., 2003, using a Daubechie's wavelet, but its robustness property are slightly less interesting). The 
asymptotic normality of such estimator is established for d G M (when the number of vanishing moments 
of the wavelet function is large enough). 

Table [H H [3] and H respectively provide the results of simulations for ARIMA(1, d, 0), ARFIMA(0, d, 0), 
ARFIMA(1, d, 1) and X^^^'i) processes for several values of d (or other parameters or 9) and N . 

Conclusions of simulations: In more than 50% of cases, the estimator d^^'' provides the smallest \/MSE 

— { / 

among the 4 semiparametric estimators. The good performances of dj^ with respect to the other estima- 
tors are particularly clear for d G [0.5, 1.20]. However, we can remark that the estimator d)^ is slightly 
less efficient for ARFIMA(1, d, 1) processes and especially X'-'^'^^ processes, that is to say, when the spectral 
density is slightly less smooth {i.e. close to a power law function). 

Comparison of MIR tests Sn and T/v with other famous stationarity or nonstationarity tests 

Monte-Carlo experiments were done for evaluating the performances of new tests Sn and T/v and for com- 
paring them to most famous stationarity (KPSS and LMC) or nonstationarity (ADF and PP) tests (see 
more details on these tests in the previous section). For these tests we selected a classical "optimal" lag 
parameter k such as: 



• fc = 



for KPSS test; 



fc = for LMC test; 
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• k = 



(n- 1)1/3 



for ADF test; 



4 (mi) 



1/4- 



for PP test; 



We also computed the results of V/S test, which is a nonparametric long-memory test (see more details 
Giraitis et at, 2007). This test is useful when d < 0.5 but it can not really distinguish the stationary case 
d € (0,0.5) to the nonstationary case d > 0.5. 

The results of these simulations with a type I error classically chosen to 0.05 are provided in Tables [SJ IH [7] 
andB 



Conclusions of simulations: For ARIMA(p, d, 0) processes with d = or d = 1 (notably AR(1) pro- 
cess when d = 0), ADF and PP tests are more accurate than our adaptive MIR tests. This is not really a 
surprise since those classical tests are unit root tests and are built from ARIMA processes. What is more 
surprising is that our MIR tests provide better results than KPSS and LMC tests (for LMC only in the case 
d= 1). 

In case of ARFIMA or X^'^'^^ processes which are stationary long-memory processes when < d < 0.5 and 
nonstationary processes when d > 0.5, all the classical tests are not able to distinguish between the cases 
d = and d = 1. Hence, if < d < 1, the asymptotic bahaviors of KPSS or LMC tests are the same 
than for d = 0: they accept the stationarity assumption Hq. In contrast, ADF and PP tests accept the 
nonstationarity assumption Hq even if the processes are stationary when < d < 0.5. This is not really a 
problem with the V/S test since this test is built for distinguishing between short and long memory processes. 
Conversely, the adaptive MIR tests (especially nonstationarity test T^) which are based on the estimation 
of the memory parameter are very efhcient for distinguishing between the stationary case d < 0.5 and the 
nonstationary case d > 0.5. 

Finally, these simulations show that the adaptive MIR tests (especially nonstationarity test T/v) are the best 
trade-off for detecting the stationarity/nonstationarity of a process. 

5.3 Application to the the Stocks and the Exchange Rate Markets 

We applied the adaptive MIR statistics as well as the other famous long-memory estimators and stationarity 
tests to Econometric data, the Stocks and Exchange Rate Markets. More precisely, the 5 following daily 
closing value time series are considered: 

1. The USA Dollar Exchange rate in Deusch-Mark, from 11/10/1983 to 08/04/2011 (7174 obs.). 

2. The USA Dow Jones Transportation Index, from 31/12/1964 to 08/04/2011 (12072 obs.). 

3. The USA Dow Jones Utilities Index, from 31/12/1964 to 08/04/2011 (12072 obs.). 

4. The USA Nasdaq Industrials Index, from 05/02/1971 to 08/04/2011 (10481 obs.). 

5. The Japan Nikkei225A Index, from 03/04/1950 to 8/04/2011 (15920 obs.). 

We considered the log-return of this data and tried to test their stationarity properties. Since station- 
arity or nonstationarity tests are not able to detect (offline) changes, we first used an algorithm de- 
veloped by M. Lavielle for detecting changes (this free software can be downloaded from his homepage: 
http : //www . math . u-psud . f r/~lavielle/prograiniiies— lavielle . htinl|) . This algorithm provides the choice 
of detecting changes in mean, in variance, and we chose to detect parametric changes in the distribution. 
Note that the number of changes is also estimated since this algorithm is based on the minimization of a 
penalized contrast. We obtained for each time series an estimated number of changes equal to 2 which are 
the following: 
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• Two breaks points for the US doUar-Deutsch Mark Exchange rate return are estimated, corresponding 
to the dates: 21/08/2006 and 24/12/2007. The Financial crisis of 2007-2011, fohowed by the late 2000s 
recession and the 2010 European sovereign debt crisis can cause such breaks. 

• Both the breaks points estimated for the US Dow Jones Transportation Index return, of the New- York 
Stock Market, correspond to the dates: 17/11/1969 and 15/09/1997. The first break change can be a 
consequence on transportation companies difficulties the American Viet-Nam war against communist 
block. The second change point can be viewed as a contagion by the spread of the Thai crisis in 1997 
to other countries and mainly the US stock Market. 

• Both the breaks points estimated for the US Dow Jones Utilities Index return correspond to the dates: 
02/06/1969 and 14/07/1998. The same arguments as above can justify the first break. The second 
at 1998 is probably a consequence of "f/ie long very acute crisis in the bond markets,..., the dramatic 
fiscal crisis and Russian Flight to quality caused by it, may have been warning the largest known by 
the global financial system: we never went too close to a definitive breakdown of relations between the 
various financial msfrt/ments "(Wikipedia). 

• The two breaks points for the US Nasdaq Industrials Index return correspond to the dates: 17/07/1998 
and 27/12/2002. The first break at 1998 is explained by the Russian flight to quality as above. The 
second break at 2002 corresponds to the Brazilian public debt crisis of 2002 toward foreign owners 
(mainly the U.S. and the IMF) which implicitly assigns a default of payment probability close to 100% 
with a direct impact on the financial markets indexes as the Nasdaq. 

• Both the breaks points estimated for the Japanese Nikkei225A Index return corresponds to the dates 
29/10/1975 and 12/02/1990, perhaps as consequence of the strong dependency of Japan to the middle 
east Oil following 1974 or anticipating 1990 oil crisis. The credit crunch which is seen as a major factor 
in the U.S. recession of 1990-91 can play a role in the second break point. 

Data and estimated instant breaks can be seen on Figure [TJ Then, we applied the estimators and tests 
described in the previous subsection on trajectories obtained in each stages for the 5 economic time series. 
These applications were done on the log-returns, their absolute values, their squared values and their 0-power 
laws with maximized for each LRD estimators. The results of these numerical experiments can be seen in 
Tables HUH 

Conclusions of numerical experiments: We exhibited again the well known result: the log-returns are 
stationary and short memory processes while absolute values or power 9 of log-returns are generally sta- 
tionary but long memory processes (for this conclusion, we essentially consider the results of Sn, Tn and 
V/S tests since the other tests have been shown not to be relevant in the cases of long-memory processes). 
However the last and third estimated stage of each time series provides generally the largest estimated values 
of the memory parameter d (for power law of log-returns) which are close to 0.5; hence, for Nasdaq time 
series, we accepted the nonstationarity assumption. 

6 Proofs 

Proof of Proposition\7\ This proposition is based on results of Surgailis et al. (2008) and was already proved 
in Bardet et Dola (2012) in the case -0.5 < d < 0.5. 

Mutatis mutandis, the case 0.5 < d < 1.25 can be treated exactly following the same steps. 

The only new proof which has to be established concerns the case d = 0.5 since Surgailis et al. (2008) do not 

provide a CLT satisfied by the (unidimensional) statistic IRN{m) in this case. Let Ym{j) the standardized 
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Daily Dow Jones Transportatio 




Daily Dow Jones Utilities Index Closings 
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Figure 1: The financial data (DowJonesTransportations, DowJonesUtilities, Nasdaqindustrials, Nikkei225A 
and US Dollar vs Deutsch Mark): original data (left) and log- return with their both estimated breaks instants 
(right) occurred at the distribution Changes 



12 



process defined Surgailis et al. (2008). Then, for d — 0.5, 



Vj>l, |7„0-)l = |E(y™(i)i;„(0))h [ cos{jx)x{c^ + 0{x^))—^dx 







sin^(f) 



Denote ^m{j) = Pm{j) = yri^i +-^2) as in (5.39) of Surgailis et al. (2008). Both inequalities (5.41) and 
(5.42) remain true for d = 0.5 and 

777"^ m"^ ^ vn 

|/i|<C— , \h\<C— =^ \i,+i,\<c— =^ \j^(j)\^\p„^(^j)\< (\i,+i,\)<c- 
3 3 3 ^ m 3 

Now let r],n{j) \Y^lfl^Yllt+^^ ^(^m(i), ^m(j + rn)^ . The Hermite rank of the function ip \s 2 and 
therefore the equation (5.23) of Surgaihs et al. (2008) obtained from Arcones Lemma remains vahd. Hence: 

9 

777 

Cov(r?„(0),r,^(j)) <C— 
from Lemma (8.2) and then the equations (5.28-5.31) remain valid for all d G [0.5, 1.25). Then for d = 0.5, 

\ -{lRN{m)-^[lRN{m)\) A AA(0, ^2(0.5)), 

with (t2(0.5) ~ (0.2524)2. □ 
Proof of Provertv \2.1\ As in Surgailis et al (2008), we can write: 

E[^^^^M]-E(^,,^-^).A(-) w.th ^:=l-2^^--^g^ 

Therefore an expansion of Rm/V^ provides an expansion of E[/i?jv(m)] when m — > oo. 

Step 1 Let / satisfy Assumption IG{d, (3) . Then we are going to establish that there exist positive real 
numbers Ci, C2 and C3 specified in (|6.1|) . (|6.2p and (|6.3p such that for 0.5 < d < 1.5 and with p{d) defined 
in (Ha, 



1. if/3<2d-l, ^ = p{d)+Ci{2-2d,P)m-'^ + 0(m~^ +m-^'^ 

2. if/3 = 2d-l, ^=p(d)+C2(2-2d,/3)m-'^+o(m-2+TO-2-/3iog(TO) + m-2'3 

3. if 2d - 1< /3 < 2d + 1, ^ = p(d) + C3(2 - 2d, P)m'l^ + o(m-^~<^ + m-^d-i j^g^^^ ^ ^j^-20A . 
4 - ^ 



if /3 = 2d + 1, ■y^= p{d) + o(m-2'i-i fog(m) + . 



Under Assumption IG{d,/3) and with Jj{a,m) defined in (j6.7p in Lemma 16.71 it is clear that, 

Rra_ _ M2 " 2d, m) + ^M2 - 2d + /?, m) + 0( J6(2 - 2d + /3 + e)) 
V;^ -1-2 ^^^2 _ 2d, m) + fi J4(2 -2d + 13, m) + 0{Ji{2 - 2d + /3 + e)) ' 

since / 0{x^-^'^+^"'^) )-^dx = 0{J,{2 - 2d + P + e)). Now using the results of Lemma O and 

Jo sin (f ) 

constants Cj£, Cj^ and C"^, j = 4, 6, £ = 1, 2 defined in Lemma 15771 



13 



1. Let < ^ < 2rf - 1< 2, i.e. -1< 2 - 2rf + /3 < 1. Then 



" "^"^ C4i(2 - 2d)rni+2''+0(rn2rf-i) +f^C4i(2 -2d + ^)TOi+2rf-/3+0(m2rf-i-/3) 

ci(74i(2-2d + /3) 



:1- 



C4i(2-2d) 



Cei (2 - 2d)+— Cei (2 2d + /S)™"'^ 

Co 



-m 



-/3 



+0(m-^) 



_^ 2C6i(2-2d) ci [C6i(2-2rf)C4i(2-2d + ^) C6i(2 - 2(i + /3) 



C41 (2 - 2d) Co L C41 (2 - 2d)C4i (2 - 2d) C41 (2 - 2d) 

As a consequence,, 

R 



Co C4i(2-2d) 

m-'5+0(m-2+m-2'3). 



-f- = p(d) + Ci(2 - 2d,^) m-^ + o(m-'^ + m.-^'^) {m 00), with < /3 < 2d - 1< 2 and 

Ci (2 - 2d, /3) := 2 ^ [^61 (2 - 2d)C4i (2 - 2d + /3) - Cgi (2 - 2d + l3)C4i (2 - 2d)] , (6.1) 

Co (^41 V-^ ~ -^"j 

and numerical experiments proves that Ci(2 — 2d,f3)/ci is negative for any d G (0.5, 1.5) and /3 > 0. 



2. Let /3 = 2d - 1, ie. 2 - 2d + /3 = 1. Then, 

_ G>i(2 - 2(/) ;>r^-2''+0(;,,2^^-i) +l^roi(l)»>i-^'^+0(l(,g(m)) 
~ C4i(2 - 2d)rni+2d+0(m2d-i) +£iC'4i(l)TOi-2d+0(log(rn)) 



= 1- 



C4i(2 - 2d) 



C6i(2-2d)+-C"6i(l)mi 

Co 



Ci C"4l(l) i_2d 
1 :::: — 7z TT-rrfn 



+0(m-2 + m"^'*-^ log(m)) 



^ 2C6i(2-2d) ci r C6i(2-2d)C'4i(l) C'^ijl) 
C4i(2 - 2d) Co l-C4i(2 - 2d)C4i(2 - 2d) C4i(2 - 2d) 



coC4i(2-2d) 

mi-2'*+0(m-2 + TO-^'^-i log(m) + m^-^-^) . 



As a consequence. 



R 



^2 _ p(d) + C2(2-2d, /3) m-''+ o(m-^+m-'^-'^ log(m)+m-^'') (m -s- 00), with < /3 = 2d - 1< 2 and 

[C6i(2-2d)C'4i(l)-C'6i(l)C4i(2-2d)], (6.2) 



(72(2-2d,/3):=2^-^ ^ 

Co <-^/n 



coC|i(2-2d) 

and numerical experiments proves that €2(2 — 2d, /3)/ci is negative for any d e [0.5, 1.5) and /3 > 0. 



3. Let 2d - 1 < /3 < 2d+ 1, i.e. 1< 2 - 2d + /3 < 3. Then, 



:l-2- 



C6i(2 - 2d)mi+2''+f C'6i(2 - 2d + ;3)mi+2'^-'5+0(mi+2'^-'^-<^ + log(m)) 



Cii{2 - 2d)TOi+2'i+fiC"4i(2 - 2d + /3)mi+2'i-/3+C)(mi+2'i-/3-<^ + m-2'i-i log(TO)) 



C4i(2 - 2d) 



C6i(2 - 2d)+— C"6i(2 - 2d + I3)m- 

Co 



^ ciC"4i(2-2d + /3) 
1 „ — ^-r — -m 



+0{m-^-'' + m-2<*-i log(m)) 



_^ 2C6i(2-2d) ^^ci rC6i(2-2d)C'4i(2-2d + /3) C'6i(2 - 2d + /3) 



Co C4i(2-2d) 

m-'^+0{m-^-'' + m-^'^-^ log(m)) . 



C41 (2 - 2d) CO L C41 (2 - 2d)C4i (2 - 2d) C41 (2 - 2d) 

As a consequence, 

^=p(d) + C3(2-2d,/3) m-f^ + o(m^'^-<^ + m-2''-Mog(m) + m-2'') (m ^ 00), 



and 



C3(2 - 2d, /3) := 2 [C6i(2 - 2d)C"4i(2 - 2d+ /3) - C'6i(2 - 2d + /3)(74i(2 - 2d)] , (6.3) 

Co L'4ll^ ^ ^"j 



and numerical experiments proves that C3(2 — 2d, /3)/ci is negative for any d G [0.5, 1.5) and /3 > 0. 
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4. Let (3 = 2d + 1. Then, Once again with Lemma 15771 



Rrr^_ _ Cei(2-2d) mi+^'^+0(m^^-i) +f^C'62(3) log(m)+0(l) 
" C4i(2 - 2d)mi+2rf+0(m2'i-i) +f-C"42(3) log(m)+0(l) 



1- 



Cii{2-2d) 



Cei{2 - 2d)+-C'62(3)m-^ log(m) 

Co 



^ 2C6i(2-2d) ^^ci r C6i(2-2d)C'42(3) 



C4i(2-2d) Co LC4i(2-2d)C4i(2-2d) C4i(2-2d) 
As a consequence, 



Ci C"42(3) _o " 

^-^ C4i(2-2d) "^ 
C'62(3) ^^-0Yog{m)+O{m-^) 



~0( 



-2d-l 



T/2 



p{d) + 0(m"^'^"^ log(m)+TO"2) (m ^ cx)), with 2 < /3 = 2(i + 1 < 4. 



(6.4) 



Step 2: A Taylor expansion of A(-) around p{d) provides: 



9A 



dp 



dp^ 



{P{d)){y^-p{d))' 



Note that numerical experiments show that 



dp 



(p) > 0.2 for any p £ (—1,1). As a consequence, using 

the previous expansions of Rm/V^ obtained in Step 1 and since E[/i?Ar(m)] = A(i?,„/y^), then for all 
< /3 < 2: 



E[Ji?Ar(m)] - Ao(d) + 



ciC'i{d,P)m-f^ + 0{m-^ + m'^f^) if /3 < 2d - 1 

ciC"2(rf,/3)m-'5 + 0(m-2 + TO-2-/3iogTO + m^2/3^ /3 = 2d - 1 

ciC"3(d,/3)m-'5 + 0(m-'5-' + 771-2^-1102^ + 771-2/3) if 2d - 1 < /? < 2(i + 1 

0(777-2^-1102777 + 777-2) if /3^l + 2d 



with C;(d, /3) 



OA 
dp 



{p{d)) Ce{2 - 2d, P) for ^ = 1, 2, 3 and Q defined in (gT]), (g^) and (lO 



□ 



Proof of TheoremU\ Using Property O if 777 ~ CiV" with C>0and(l + 2^)-i<a<l then 
^N/m (E[/i?Ar(77i)] — Ao(d)) — ^ and it implies that the multidimensional CLT (|2.5I) can be replaced 

by 



N 

777 V 



(iRNirrij) - Md) 



l<j<p N^oo 



AA(o,rp(d)). 



(6.5) 



It remains to apply the Delta-method with the function A^^ to CLT (|6.5I) . This is possible since the 
function d — >■ Ao{d) is an increasing function such that AQ(d) > and (Ag ^)'(Ao(d)) = 1/Ag(d) > for aU 
de (—0.5, 1.5). It achieves the proof of Theorem [TJ □ 



Proof of Proposition\^ See Bardet and Dola (2012). 
Proof of Theorem\M See Bardet and Dola (2012). 



□ 
□ 



Appendix 

We first recall usual equalities frequently used in the sequel: 

Lemma 6.1. For all X > 

2 f°° sm{Xx) , 4a f°° sm^(Xx) , 
J. For a € 0,2 , — — -/ ^ — '-dx ^ —— ^—r^dx = 



|A|°-i 7o 2'»|A|° 7o x"^^ r(a)sin(if)' 

^ ^ ^ , , 1 [^siu\Xx)^ 16 sin6(Aa;)^ 2^-^\X\^-^ n 

2. For b£ {-1,1),— r— _^dx = — — — — — -^x / —^dx = 



21-" -Uo 2;4-& -15 + 6 -23-" -33-" x^-" ' 4 r(4 - 6) sin(ii 
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3. Forb G (1,3), 



1-21-'' 



siii^(Ax) 



dx 
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15 - 6 • 23-'' + a^-f- ■ 



23-fc|;^|3-fc 



dx 



4 r(4-6) sin(MZL)- 



Proof. These equations are given or deduced (using decompositions of sin^(-) and integration by parts) from 
(see Doukhan et ai, p. 31). 



Lemma 6.2. For j ~ 4, 6, denote 



Jj{a, m) 



^^sin^ 
sin^d) 



dx. 



□ 



(6.6) 



Then, we have the following expansions when m — > cx); 



Jj(a, m) 



Cji{a)rrr-^'' + 0{m^-'') if -l< a < I 

qi(l)TO3-Q + o(iog(TO)) ifa^l 

C'jj^{a)m^^'' + 0{l) i/l<a<3 

C;.2(3) log(m) + 0(l) ifa = 3 

C"ji(a) + 0(TO-«'^-3)A2)) ^fa>3 



(6.7) 



with the following real constants (which do not vanish for any a on the corresponding set): 
4 7r(l-^) ^ _ 7r(15 - 6 • 23-'' +33-°) 



• C4i(a) := 



(3-a)r(3-a) sin(^^^) 



and CGi{a) 



3 — a 



-'{l<a<3} + 16 



4(3-a)r(3-a) sin(^5^) 
1 



and Cgi(a) := 



16 



isin^ 

'7/4— a 

y 



dy- 



1 sin^(^) 

-A^dy + 2 1 ( - 4 cos(y) + cos(2y) j dy 

y J 1 y ^ ' 

5 1 

-'{l<a<3} + T. 



3 — a 



1 



2 7i y4-a 

C42(a) := (s • i{a=3} + -^{0=1}) a"rf •^62(0) (s • J{q=3} + - • 



15 cos(j/) + 6 cos(2y) — cos(3y) ) dy 



8 Jo sm* f ) 16^0 



sin^(f) 



-dx. 



Proof. The proof of these expansions follows the steps than those of Lemma 5.1 in Bardet and Dola (2012). 
Hence we write for j = 4, 6, 



Jjia, m) 



,mx , 1 



Jj{a,m)+ I x°-sm\ — )j—^dx+ I x" sm\ — ) 



2 '(I) 



, mx , 2 1 



2 '3(1)2 



dx 



with 



Jj(a, m) 



The expansions when m — > 00 of both the right hand sided integrals in (j6.8p are obtained from Lemma |6. II 
It remains to obtain the expansion of Jj{a,m). Then, using classical trigonometric and Taylor expansions: 



sin''(|) = -^(3-4003(2/) +cos(2y)) and 



1 



1 



8 

sin6(|) = l(i0-15cos(y) + 6cos(2y)-cos(3y)) and 1 + 11 - '^"'^^^ 



12 1 11 , 



sin*(y) r 3y2 45 



3y3 sin^(2/) 945 



the expansions of Jj (a, m) can be obtained. 

Numerical experiments show that C'^i{a) 7^ 0, Cgi(a) ^ 0, (^42(0) ^ and C'^2{a) ^ 0. 



□ 
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N = 200 
ARIMA(l,d, 0) 


d = Q 
0=-O.7 


d = 

<p=-0.8 


d = 

(f>=-0.9 


d = 

(}i=-0.95 


d= 1 

0=-l 


d=l 
4>=-1.05 


d=l 
(f)=-l.l 


d= 1 

<f)=-1.2 


mean a]^ 
Vmse d^j^^^ 


0.375 
0.399 


0.518 
0.536 


0.717 
0.729 


0.822 
0.833 


0.964 
0.135* 


0.976 
0.129* 


0.980 
0.134* 


1.008 
0.127* 


mean dj^ ' 

Vmse df-^^^ 


0.244 
0.286* 


0.426 
0.451 


0.673 
0.689 


0.858 
0.875 


1.057 
0.173 


1.075 
0.178 


1.069 
0.179 


1.083 
0.189 


mean a]^ 


0.003 

0.361 


0.124 
0.380* 


0.312 
0.519* 


0.623 
0.739* 


0.946 
0.344 


0.934 
0.349 


0.977 
0.314 


0.947 
0.332 


mean a]y 


0.461 
0.491 


0.600 
0.615 


0.747 
0.756 


0.831 
0.843 


0.902 
0.175 


0.925 
0.163 


0.933 
0.162 


0.968 
0.159 



N = 500 
ARIMA(1, d, 0) 


d = Q 
0=-O.7 


d = 
0=-O.8 


d = 
0=-O.9 


d = 
0=-O.95 


d= 1 
0=-l 


d = 1 
</>=-1.05 


d= 1 
(j>=-l.l 


d= 1 
(j>=-l.2 


mean aj^ 

Vmse c?^^' 


0.232 

0.279 


0.355 

0.400 


0.629 

0.649 


0.804 

0.811 


0.980 
0.103* 


0.981 
0.101* 


0.976 
0.094* 


0.995 
0.104* 


mean d j^ 
VMSEd^^^^ 


0.068 
0.154* 


0.157 
0.212 


0.390 
0.412 


0.644 
0.659 


1.057 
0.174 


1.057 
0.164 


1.055 
0.171 


1.060 
0.174 


mean d]^ 

Vmse ^f^^' 


0.044 
0.176 


0.092 
0.209* 


0.275 
0.321* 


0.560 
0.589* 


0.976 
0.164 


0.969 
0.159 


0.976 
0.169 


0.973 
0.163 


■:iWAV) 

mean d j^ 

VmseS-^^^^ 


0.291 
0.337 


0.442 
0.479 


0.679 
0.694 


0.818 
0.825 


0.940 
0.118 


0.944 
0.116 


0.948 
0.108 


0.979 
0.106 



N = 1000 
ARIMA(l,d, 0) 


d = 
(p=-0.7 


d = 
(p=-0.8 


d = 
<^=-0.9 


d = 
(p=-0.95 


d = 1 
0=-l 


d = 1 
(^=-1.05 


d = 1 
(f)=-l.l 


d = 1 
(f)=-1.2 


mean dj^ 

Vmse <?^-^' 


0.171 
0.221 


0.284 
0.338 


0.579 
0.607 


0.774 
0.784 


0.975 
0.074* 


0.981 
0.073* 


0.994 
0.067* 


0.995 
0.073* 


liMS) 

mean d j^ 

Vmse d</^' 


0.057 
0.108* 


0.137 

0.162 


0.340 
0.351 


0.574 

0.580 


1.027 

0.103 


1.031 

0.103 


1.029 

0.102 


1.031 

0.099 


mean d]^ 
x/MSEd<^°^' 


0.035 
0.112 


0.091 
0.146* 


0.266 
0.289* 


0.515 
0.528* 


0.987 
0.096 


0.994 
0.100 


0.998 
0.102 


0.987 
0.095 


mean dj^ 
7M5Ed<7^^) 


0.165 
0.237 


0.252 
0.306 


0.424 
0.467 


0.663 
0.689 


0.949 
0.090 


0.958 
0.092 


0.963 
0.097 


0.976 
0.090 



N = 5000 
ARFIMA(l,d,0) 


d = 
(f>=-0.7 


d = 
(p=-0.8 


d = 
<^=-0.9 


d = 
(p=-0.95 


d = 1 
0=-l 


d = 1 
(^=-1.05 


d = 1 
<p=-l.l 


d = 1 
<^i=-1.2 


mean dj^ 

Vmse 4^^' 


0.064 
0.100 


0.083 
0.131 


0.212 
0.303 


0.661 
0.714 


0.989 
0.034* 


0.989 
0.035* 


0.996 
0.032* 


1.003 
0.036 


mean d]y 

Vmse d</^' 


0.014 
0.042* 


0.061 
0.075* 


0.219 
0.224 


0.423 
0.426 


1.006 
0.044 


1.002 
0.041 


1.007 
0.045 


1.005 
0.041 


mean d]y ' 

Vmse d<«°^) 


0.034 
0.053 


0.080 
0.090 


0.243 
0.247 


0.479 
0.481 


0.997 
0.038 


0.997 
0.036 


0.997 
0.037 


0.999 
0.034* 


mean d]^ 


0.031 
0.136 


0.045 
0.134 


0.047 
0.151* 


0.167 
0.212* 


0.980 
0.083 


0.982 
0.071 


0.977 
0.099 


1 

0.067 



Table 1: : Comparison between d^^' and other famous semiparametric estimators of d {d^^^\ and d^'*^^) 

applied to ARIMA(l,d, 0) process {Xt + <pXt-i = et), with several <f) and N values 
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JV = 200 
ARFIMA(0,d,0) 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


d = OA 


d = 0.5 


Od= .6 


d = 0.8 


d=l 


d= 1.2 


mecui 


-0.267 
0.164 


-0.134 
0.131 


0.024 
0.115 


0.183 
0.135 


0.365 
0.122* 


0.467 
0.115* 


0.567 
0.129* 


0.770 
0.124* 


0.971 
0.139* 


1.179 
0.133* 


mean a]^ ' 


-0.354 
0.166 


-0.170 
0.151 


0.008 
0.144 


0.201 
0.163 


0.403 
0.152 


0.536 
0.151 


0.644 
0.168 


0.883 
0.192 


1.078 
0.188 


1.181 
0.167 


mean dj^ 

Vmse 


-0.310 
0.261 


-0.159 
0.301 


-0.012 
0.361 


0.201 
0.382 


0.357 
0.411 


0.476 
0.368 


0.604 
0.374 


0.781 
0.369 


0.984 
0.361 


1.065 
0.319 


mean aj^ 


-0.408 
0.094* 


-0.229 
0.107* 


-0.040 
0.101* 


0.143 
0.120 


0.352 
0.124 


0.447 
0.135 


0.531 
0.144 


0.715 
0.160 


0.915 
0.167 


1.091 
0.201 



N = 500 
ARFIMA(0,d,0) 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


d = .6 


d = 0.8 


d = 1 


d = 1.2 


mean d j^ 


-0.287 

0.135 


-0.163 

0.089 


-0.007 

0.092 


0.187 

0.096 


0.379 
0.094* 


0.464 
0.100* 


0.558 
0.100* 


0.778 
0.100* 


0.963 
0.111* 


1.180 
0.114* 


mean a]^ 
^MSEd^^^^ 


-0.344 
0.157 


-0.181 
0.144 


0.001 
0.135 


0.209 
0.141 


0.411 
0.158 


0.522 
0.154 


0.636 
0.168 


0.870 
0.183 


1.034 
0.173 


1.158 
0.172 


-iROB) 

incaii 


-0.344 
0.146 


-0.206 
0.174 


-0.029 
0.164 


0.196 
0.172 


0.376 
0.190 


0.486 
0.182 


0.596 
0.187 


0.799 
0.176 


0.956 
0.168 


1.058 
0.214 


-XiWAV) 

mean a]^ 


-0.403 
0.065* 


-0.218 
0.063* 


-0.027 
0.067* 


0.166 
0.080* 


0.356 
0.107 


0.447 
0.107 


0.537 
0.108 


0.753 
0.109 


0.928 
0.124 


1.134 
0.137 



N = 1000 
ARFIMA(0,d,0) 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


Od = .6 


d = 0.8 


d = 1 


d = 1.2 


mean a]^ 
VMSE ^^^^ 


-0.311 
0.105 


-0.170 
0.065 


-0.001 
0.063 


0.181 
0.059* 


0.374 
0.065* 


0.476 
0.068* 


0.575 
0.070* 


0.774 
0.075* 


0.973 
0.074* 


1.181 
0.079* 


mean dS^^^ 

Vmse d^/^' 


-0.360 

0.106 


-0.183 

0.092 


-0.004 

0.086 


0.200 

0.090 


0.398 

0.090 


0.506 

0.092 


0.628 
0.108 


0.843 

0.109 


1.018 

0.106 


1.116 

0.129 


mean a]^ 

Vmse 


-0.380 
0.097 


-0.204 
0.115 


-0.012 
0.112 


0.192 
0.106 


0.388 
0.113 


0.501 
0.107 


0.617 
0.123 


0.811 
0.105 


0.982 
0.104 


1.075 
0.156 


mean a]^ 


-0.407 
0.059* 


-0.213 
0.050* 


-0.026 
0.052* 


0.164 
0.067 


0.357 
0.090 


0.460 
0.088 


0.550 
0.095 


0.761 
0.089 


0.938 
0.114 


1.137 
0.103 



N = 5000 
ARFIMA(0,d,0) 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


Od = .6 


d = 0.8 


d = 1 


d = 1.2 


mean dj^ 

Vmse 4^^' 


-0.339 
0.070 


-0.177 
0.040* 


-0.003 
0.031* 


0.189 
0.033* 


0.386 
0.035* 


0.483 
0.035* 


0.589 
0.035* 


0.786 
0.035* 


0.984 
0.042 


1.199 
0.052* 


mean d]y 

Vmse d</^> 


-0.385 
0.043 


-0.195 
0.044 


0.001 
0.044 


0.197 
0.042 


0.400 
0.043 


0.502 
0.044 


0.612 
0.051 


0.825 
0.059 


1.006 
0.043 


1.078 
0.142 


mean dj^ 


-0.395 
0.042* 


-0.203 
0.043 


-0.002 
0.042 


0.197 
0.041 


0.402 
0.041 


0.504 
0.040 


0.609 
0.045 


0.820 
0.051 


0.997 
0.035* 


1.070 
0.147 


yfWAV) 
mean a]y 


-0.416 
0.080 


-0.210 
0.055 


-0.025 
0.083 


0.175 
0.078 


0.372 
0.082 


0.476 
0.074 


0.575 
0.086 


0.782 
0.068 


0.979 
0.083 


1.171 
0.103 



Table 2: : Comparison between c% and other famous semiparametric estimators of d (d}^ \ '^n and d]y, ) 
applied to ARFIMA(0, d, 0) process, with several d and N values 
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N = 200 
ARFIMA(l,d,l) 
(/>= -0.3 ; 6» = 0.7 


d= -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


Od= .6 


d = 0.8 


d = l 


d = 1.2 


mean a., ' 

N 

yjMSE 4^-"' 


-0.154 
0.273 


-0.009 

0.231 


0.151 

0.197 


0.310 
0.176* 


0.476 
0.150* 


0.576 
0.147* 


0.669 
0.146* 


0.861 
0.137* 


1.051 
0.143* 


1.234 
0.129* 


mean d]^ 

^Tmse 


-0.241 
0.223* 


-0.055 
0.208* 


0.124 
0.190* 


0.313 
0.200 


0.503 
0.183 


0.627 
0.197 


0.719 
0.202 


0.925 
0.215 


1.104 
0.199 


1.196 
0.169 


mean a]y 


-0.312 
0.259 


-0.152 
0.304 


-0.006 
0.361 


0.210 
0.383 


0.364 
0.411 


0.483 
0.368 


0.610 
0.375 


0.787 
0.368 


0.988 
0.360 


1.068 
0.317 


mean a]^ 


-0.199 
0.253 


-0.012 
0.236 


0.182 
0.232 


0.364 
0.232 


0.561 
0.214 


0.635 
0.203 


0.731 
0.195 


0.890 
0.184 


1.039 
0.157 


1.183 
0.168 



N = 500 

ARFIMA(l,d,l) 
<j) = -0.3 ; 6* = 0.7 


d= -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


0d= .6 


ci = 0.8 


d=l 


d = 1.2 


mean a]^ 
^/MSE 4^^^ 


-0.216 

0.204 


-0.090 

0.147 


0.062 
0.122* 


0.255 
0.128* 


0.435 
0.118* 


0.529 
0.116* 


0.617 
0.114* 


0.827 
0.114* 


0.999 

0.1161 


1.214 
0.113* 


mean a]y 

Vmse 


-0.363 
0.149 


-0.197 
0.144* 


-0.015 
0.135 


0.194 
0.141 


0.398 
0.158 


0.511 
0.153 


0.626 
0.165 


0.866 
0.182 


1.036 
0.172 


1.160 
0.171 


mean d\, 
VMSEd^^°^'^ 


-0.346 
0.1450* 


-0.203 
0.175 


-0.022 
0.162 


0.202 
0.172 


0.382 
0.189 


0.492 
0.182 


0.602 
0.187 


0.804 
0.177 


0.959 
0.167 


1.060 
0.212 


mean rf}^ 


-0.292 
0.170 


-0.090 
0.169 


0.093 
0.157 


0.303 
0.158 


0.502 
0.164 


0.579 
0.157 


0.679 
0.151 


0.864 
0.142 


1.017 
0.112* 


1.190 
0.121 



= 1000 
ARFIMA{l,d,l) 

<f> = -0.3 ; 6» = 0.7 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


0d= .6 


ci = 0.8 


d=l 


d = 1.2 


mean aj^ 
^/MSE 


-0.245 
0.166 


-0.110 
0.121 


0.055 
0.104 


0.237 
0.095 


0.425 
0.091* 


0.528 
0.091* 


0.628 
0.092* 


0.815 
0.093* 


1.005 
0.090* 


1.212 
0.086* 


i(mS) 
mean a]y 

^/MSE 4^^) 


-0.376 
0.099 


-0.197 
0.090* 


-0.018 
0.088* 


0.186 
0.091* 


0.386 
0.092 


0.495 
0.093 


0.619 
0.108 


0.840 
0.108 


1.017 
0.105 


1.117 
0.129 


mean a]^ 


-0.379 
0.096* 


-0.199 

0.116 


-0.007 

0.111 


0.197 

0.105 


0.393 
0.112 


0.507 

0.107 


0.622 

0.124 


0.815 
0.105 


0.986 

0.103 


1.076 

0.155 


JCWAV) 

mean a]^ 


-0.334 
0.134 


-0.136 
0.129 


0.046 
0.129 


0.247 
0.126 


0.459 
0.132 


0.555 
0.137 


0.645 
0.128 


0.837 
0.126 


1.006 
0.103 


1.183 
0.099 



N = 5000 
ARFIMA(l,d,l) 
(t>= -0.3 ; 6» = 0.7 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


0d= .6 


d = 0.8 


d=l 


d = 1.2 


JylR) 

mean aj^ 
^/MSE 4^^' 


-0.298 
0.111 


-0.143 
0.075 


0.037 
0.064 


0.226 
0.062 


0.423 
0.058 


0.523 
0.051 


0.621 
0.051 


0.813 
0.048* 


1.005 
0.054 


1.219 
0.063* 


mean d]^ ' 

Vmse 4^^^' 


-0.391 
0.042* 


-0.200 

0.044 


-0.004 

0.044 


0.192 

0.043 


0.396 

0.043 


0.498 

0.044 


0.609 

0.051 


0.824 

0.059 


1.006 

0.043 


1.078 

0.142 


%ROB) 
mean dj^ 

VMSEd^^"^^ 


-0.392 
0.042 


-0.198 
0.043* 


0.003 
0.042* 


0.201 
0.041* 


0.407 
0.042* 


0.509 
0.041* 


0.613 
0.046* 


0.825 
0.053 


1.002 
0.035* 


1.072 
0.146 


liWAV) 
mean dff 

VMSEd^^^^^ 


-0.396 
0.085 


-0.187 
0.094 


0.003 
0.097 


0.196 
0.107 


0.405 
0.110 


0.519 
0.100 


0.605 
0.124 


0.818 
0.102 


1.014 
0.107 


1.192 
0.101 



Table 3: : Comparison between c% and other famous semiparametric estimators of d {d}^ i ^jv ^^'^ 'vv ) 
applied to ARFIMA(l,d, 1) process (with <p = —0.3 and 9 = 0.7), with several d and N values. 
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N = 200 


d = -0.4 


d= -0.2 


d = 


d = 0.2 


d = OA 


d = 0.5 


Od= .6 


d = 0.8 


d= 1 


d= 1.2 


mean a]^^ ' 
\/MSE S^^' 


-0.328 
0.104* 


-0.221 
0.107* 


-0.097 
0.155* 


0.080 
0.173* 


0.250 
0.194 


0.320 
0.219 


0.435 
0.212 


0.623 
0.223 


0.852 
0.194 


1.021 
0.224 


mean a]^ 


-0.455 
0.165 


-0.266 - 
0.172 


0.080 
0.173 


0.122 
0.173 


0.320 
0.172* 


0.409 
0.183* 


0.527 
0.178* 


0.733 
0.191* 


0.990 
0.185* 


1.219 
0.183* 


mean a]^ 


-0.322 
0.251 


-0.236 
0.277 


-0.054 
0.339 


0.118 
0.362 


0.330 
0.393 


0.411 
0.381 


0.494 
0.396 


0.727 
0.376 


0.944 
0.402 


1.120 
0.328 


mean dj^ ' 


-0.636 
0.259 


-0.424 
0.254 


-0.228 
0.254 


-0.014 
0.246 


0.184 
0.244 


0.253 
0.283 


0.358 
0.279 


0.552 
0.285 


0.780 
0.266 


0.992 
0.251 



N = 500 


d=-OA 


d = -0.2 


d = 


d = 0.2 


d = OA 


d = 0.5 


Od = .6 


d = 0.8 


d=l 


d = 1.2 


mean a]^ 
Vmse <?^^' 


-0.358 
0.074* 


-0.257 
0.092* 


-0.096 
0.129* 


0.087 
0.153 


0.275 
0.161 


0.376 
0.161 


0.470 
0.167 


0.677 
0.167* 


0.888 
0.154* 


1.079 
0.162* 


mean a]^ 
VMSEd^^^^ 


-0.373 
0.163 


-0.221 
0.134 


-0.045 
0.150 


0.167 
0.131* 


0.363 
0.148* 


0.480 
0.146* 


0.571 
0.157* 


0.789 
0.175 


1.023 
0.173 


1.265 
0.204 


mean d)^ 

Vmse 


-0.408 
0.113 


-0.247 
0.162 


-0.071 
0.185 


0.170 
0.159 


0.351 
0.178 


0.454 
0.188 


0.547 
0.184 


0.742 
0.189 


0.9530 
0.182 


1.1573 
0.188 


J^WAV) 

mean a]y 


-0.593 
0.214 


-0.387 
0.207 


-0.188 
0.208 


0.030 
0.196 


0.233 
0.199 


0.339 
0.196 


0.430 
0.203 


0.633 
0.196 


0.8380 
0.199 


1.0589 
0.183 



N = 1000 


d= -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


Od = .6 


d = 0.8 


d= 1 


d = 1.2 


mean d)^ ' 
y/MSE <?j^^' 


-0.382 
0.044* 


-0.253 
0.081* 


-0.096 
0.120 


0.081 
0.138 


0.296 
0.127 


0.376 
0.148 


0.485 
0.137 


0.686 
0.140 


0.8912 
0.134 


1.0898 
0.132 


mean a]^ 
^MSE 4^^' 


-0.409 

0.092 


-0.221 

0.093 


-0.035 
0.101* 


0.149 
0.107* 


0.361 
0.103* 


0.466 
0.102* 


0.565 
0.103* 


0.775 
0.111* 


0.996 
0.105* 


1.219 

0.122 


mean a]^ 


-0.420 
0.077 


-0.244 
0.110 


-0.037 
0.110 


0.150 
0.125 


0.358 
0.119 


0.459 
0.117 


0.562 
0.120 


0.754 
0.128 


0.959 
0.112 


1.164 
0.119* 


mean a]^ ' 


-0.564 
0.182 


-0.350 
0.173 


-0.152 
0.176 


0.051 
0.172 


0.263 
0.164 


0.363 
0.169 


0.478 
0.160 


0.674 
0.163 


0.875 
0.161 


1.089 
0.150 



N = 5000 


d= -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


Od = .6 


d = 0.8 


d = 1 


d = 1.2 


KIR) 
mean d\r 

Vmse d'^"" 


-0..392 
0.028' 


-0.257 
().0(),S 


-0.089 
0.100 


0.110 
0.102 


0.324 
0.091 


0.417 
0.099 


0.533 
0.087 


0.730 
0.08.-) 


0.939 
0.078 


1.138 
0.080 


TIMS) 
mean a]^ ' 

Vmse 


-0.409 
0.047 


-0.216 
0.046* 


-0.022 
0.048* 


0.175 
0.051* 


0.374 
0.051* 


0.473 
0.056 


0.577 
0.050 


0.779 
0.052 


0.983 
0.048* 


1.190 
0.052 


mean d]^ 


-0.426 

0.048 


-0.225 

0.049 


-0.029 

0.052 


0.172 

0.054 


0.372 

0.052 


0.469 
0.056* 


0.574 
0.050* 


0.773 
0.049* 


0.969 

0.053 


1.173 
0.051* 


mean d]y 


-0.493 
0.116 


-0.293 
0.118 


-0.093 
0.110 


0.111 
0.116 


0.314 
0.118 


0.419 
0.117 


0.512 
0.137 


0.723 
0.105 


0.931 
0.102 


1.133 
0.109 



Table 4: : Comparison between and other famous semiparametric estimators of d (d}^ i and d],, ) 

applied to X^^'^^ process with several d and N values. 
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N = 200 


d = 


d = 


d = 


d = 


d=l 


d=l 


d=l 


d=l 


ARIMA(l,d,0) 


(p=-0.7 


(p=-O.S 


0=-O.9 


0=-O.95 


<P=-1 


0=-l.O5 


</>=-l.l 


4>=-1.2 


Sn- Accepted Ho 


0.997 


0.870 


0.303 


0.130 


0.020 


0.027 


0.023 


0.010 


KPSS: Accepted Hq 


0.637 


0.470 


0.183 


0.087 


0.030 


0.017 


0.013 


0.017 


LMC; Accepted Ho 


0.983 


0.917 


0.750 


0.657 


0.530 


0.277 


0.077 





Tm- Rejected H'^ 


0.350 


0.070 


0.003 

















ADF: Rejected if^ 


1 


1 


0.957 


0.677 


0.047 


0.037 


0.063 


0.050 


PP : Rejected J?^ 


1 


1 


1 


0.797 


0.005 


0.047 


0.060 


0.037 


V/S : Accepted Ho 


0.7-57 


0.527 


0.2.33 


0.090 


0.017 


0.017 


0.007 


0.013 



N = 500 


d = 


d = 


d = 


d = 


d=l 


d=l 


d=l 


d=l 


ARIMA(l,d,0) 


<t>=-0.7 


4>=-o.s 


4>=-0.9 


0=-O.95 


4,=-l 


4>=-im 


<j)=-l.l 


<j)=-1.2 


Sn- Accepted Hq 


1 


0.950 


0.353 


0.087 


0.010 


0.007 


0.003 


0.007 


KPSS: Accepted Ho 


0.793 


0.673 


0.303 


0.117 


0.003 





0.003 


0.003 


LMC: Accepted Ho 


1 


0.993 


0.860 


0.707 


0.470 


0.140 


0.010 





Tn- Rejected H'^ 


0.787 


0.450 


0.060 

















ADF: Rejected 


1 


1 


1 


0.997 


0.067 


0.050 


0.033 


0.033 


PP : Rejected H'^^ 


1 


1 


1 


1 


0.070 


0.050 


0.030 


0.030 


V/S : Accepted Ho 


0.773 


0.653 


0.260 


0.030 





0.003 









N = 1000 


d = 


d = 


d = 


d = 


d=l 


d=l 


d=l 


d=l 


ARIMA(l,d,0) 


0=-O.7 


(l)=-0.8 


(^=-0.9 


0=-O.95 


4>=-l 


(7i=-1.05 


<j}=-l.l 


4,=-l.2 


Sf^r: Acceirted Ho 


1 


0.950 


0.377 


0.117 














KPSS: Accepted Ho 


0.817 


0.723 


0.447 


0.107 














LMC: Accepted Ho 


1 


1 


0.963 


0.797 


0.560 


0.060 








Tjv: Rejected H'^ 


0.957 


0.680 


0.173 


0.007 














ADF: Rejected H^ 


1 


1 


1 


1 


0.057 


0.047 


0.043 


0.047 


PP : Rejected H'^ 


1 


1 


1 


1 


0.057 


0.047 


0.027 


0.033 


V/S : Accepted Ho 


0.817 


0.693 


0.267 


0.037 















N = 5000 


d = 


d = 


d = 


d = 


d= 1 


d=l 


d=l 


d=l 


ARIMA(l,d,0) 


4)=-Q.7 


4>=-0.8 


4>=-0.9 


4>=-0.95 


0=-l 


4>=-im 


(j)=-l.l 


<t>=-\.2 


Sn- Accepted Ho 


1 


0.997 


0.877 


0.293 














KPSS: Accepted Ho 


0.930 


0.867 


0.687 


0.430 














LMC: Accepted Ho 


1 


1 


1 


0.963 


0.570 











Tat: Rejected 


1 


0.997 


0.850 


0.240 














ADF: Rejected H^ 


1 


1 


1 


1 


0.033 


0.057 


0.057 


0.043 


PP : Rejected H'^ 


1 


1 


1 


1 


0.037 


0.060 


0.060 


0.040 


V/S : Accepted Ho 


0.903 


0.893 


0.747 


0.183 















Tabic 5: Comparisons of stationarity and nonstationarity tests from 300 independent replications of ARIMA(1, d, 0) 
processes {Xt + 4iXt-\_ = £t) for several values of 4> and N. The accuracy of tests is measured by the frequencies of 
trajectories "accepted as stationary" (accepted Ho or rejected Ho) among the 300 replications which should be close 
to 1 for d € (-0.5, 0.5) and close to for d € [0.5, 1.2] 
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N = 200 


d = -0.4 


d= -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 


ARFIMA(0,d,0) 






















Sn- Accepted Ho 


1 


1 


1 


0.997 


0.997 


0.980 


0.803 


0.203 


0.010 


0.003 


KPSS: Accepted Ho 


1 


1 


0.950 


0.727 


0.340 


0.240 


0.097 


0.057 


0.017 


0.013 


LMC: Accepted Ho 


1 


1 


0.933 


0.360 


0.017 


0.007 





0.010 


0.483 


1 


Tn- Rejected H'^ 


1 


1 


0.997 


0.860 


0.380 


0.100 


0.033 











ADF: Rejected H'^ 


1 


1 


1 


1 


0.873 


0.553 


0.387 


0.170 


0.043 


0.010 


PP : Rejected H'^^ 


1 


1 


1 


1 


1 


0.927 


0.747 


0.360 


0.063 





V/S : Accepted Ho 


1 


1 


0.953 


0.617 


0.347 


0.173 


0.107 


0.033 


0.017 


0.007 



N = 500 


d = -0.4 


d= -0.2 


d = 


ci = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 


ARFIMA(0,d,0) 






















Sm- Accepted -ffo 


1 


1 


1 


1 


1 


0.960 


0.793 


0.070 


0.010 





KPSS: Accepted Ho 


1 


1 


0.953 


0.610 


0.167 


0.110 


0.030 


0.010 


0.003 





LMC: Accepted Ho 


1 


1 


0.943 


0.223 














0.493 


1 


Tjv: Rejected Hl^ 


1 


1 


0.997 


0.950 


0.513 


0.147 


0.023 











ADF: Rejected H[^ 


1 


1 


1 


1 


0.973 


0.737 


0.607 


0.260 


0.053 


0.010 


PP : Rejected H^^ 


1 


1 


1 


1 


1 


0.990 


0.910 


0.437 


0.053 


0.007 


V/S : Accepted Ho 


1 


1 


0.943 


0.480 


0.170 


0.053 


0.033 












N = 1000 


d = -0.4 


d= -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 


ARFIMA(0,d,0) 






















Sr>f: Accepted Ho 


1 




1 


1 


1 


0.960 


0.553 


0.020 








KPSS: Accepted Ho 


1 




0.957 


0.503 


0.100 


0.060 


0.010 


0.003 








LMC: Accepted Ho 


1 




0.943 


0.107 














0.513 


1 


T^,: Rejected H^ 


1 




1 


1 


0.773 


0.180 


0.010 











ADF: Rejected H'^ 


1 




1 


1 


0.997 


0.867 


0.747 


0.343 


0.047 


0.003 


PP : Rejected H^ 


1 




1 


1 


1 


0.997 


0.973 


0.527 


0.043 





V/S : Accepted Ho 


1 




0.960 


0.403 


0.057 


0.017 















N = 5000 


d= -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d = 1.0 


d = 1.2 


ARFIMA(0,d,0) 






















^jv^ Accepted Ho 


1 


1 


1 


1 


1 


0.980 


0.067 











KPSS: Accepted Ho 


1 


1 


0.947 


0.380 


0.010 


0.010 














LMC: Accepted Ho 


1 


1 


0.927 


0.017 














0.507 


1 


Tjv: Rejected H/, 


1 


1 


1 


1 


0.960 


0.240 














ADF: Rojcctod H,', 


1 


1 


1 


1 


1 


0.977 


0.923 


0.417 


0.067 





PP : i-tcjtxLcd Yi,', 


1 


1 


1 


1 


1 


1 


1 


O.GIT 


0.0()() 





V/S : Accepted Ho 


1 


1 


0.957 


0.153 





















Tabic 6: Comparisons of stationarity and nonstationarity tests from 300 independent replications of ARFIMA(0, d, 0) 
processes for several values of d and A''. The accuracy of tests is measured by the frequencies of trajectories "accepted 
as stationary" (accepted Ho or rejected Ho) among the 300 replications which should be close to 1 for d e (—0.5, 0.5) 
and close to for d € [0.5, 1.2] 
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TV = 200 


























W — 9 
ft — — u.z 


(I — u 


W — 9 
U — U.Z 


W — OA 


Or U.O 


a, — U.D 


Or U.O 


/7 — 10 


/7 — 19 


4> = -0.3 -,9 = 0.7 






















Sn- Accepted Ho 


1 


1 


1 


0.997 


0.930 


0.787 


0.493 


0.073 








KPSS: Accepted Ho 


1 


0.990 


0.927 


0.547 


0.317 


0.203 


0.1,37 


0.040 


0.003 


0.003 


LMC; Acc(-i)te<l /io 


0.1) l:-! 


0.810 


0.900 


1 


1 


1 


1 


1 


1 


1 


Tj\r: Rejected H^ 


1 


0.997 


0.920 


0.497 


0.100 


0.030 


0.010 











ADF: Rejected H'^ 


1 


1 


1 


1 


0.797 


0.610 


0.417 


0.203 


0.053 


0.027 


PP : Rejected i?^ 


1 


1 


1 


1 


0.933 


0.777 


0.537 


0.183 


0.027 


0.003 


V/S : Accepted Hq 


1 


0.997 


0.910 


0.537 


0.310 


0.157 


0.097 


0.033 


0.020 


0.003 



N = 500 






















ARFIMA(l,d,l) 


d = -0.4 


d = -0.2 


d = 


ci = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d = 1.0 


d = 1.2 


4> = -0.3 ; 6» = 0.7 






















Sm- Accepted Ho 


1 


1 


1 


1 


0.970 


0.793 


0.563 


0.080 


0.010 





KPSS: Accepted Hq 


1 


1 


0.913 


0.557 


0.193 


0.087 


0.060 


0.003 


0.007 





LMC: Accepted Hq 


0.967 


0.887 


1 


1 


1 


1 


1 


1 


1 


1 


Tjv: Rejected Hl^ 


1 


1 


1 


0.850 


0.207 


0.043 


0.003 











ADF: Rejected H[^ 


1 


1 


1 


1 


0.977 


0.760 


0.573 


0.203 


0.057 


0.027 


PP : Rejected H^^ 


1 


1 


1 


1 


0.997 


0.880 


0.757 


0.250 


0.013 





V/S : Accepted Hq 


1 


1 


0.917 


0.460 


0.163 


0.083 


0.030 


0.003 


0.003 






AT = 1000 






















ARFIMA(l,d,l) 


d = -0.4 


d= -0.2 


d = 


ci = 0.2 


ci = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 


<!> = -0.3 -,9 = 0.7 






















Sn- Accepted Hq 


1 


1 


1 


1 


0.980 


0.767 


0.377 


0.033 








KPSS: Accepted Ho 


1 


1 


0.937 


0.477 


0.117 


0.040 


0.010 


0.003 








LMC: Accepted Hq 


0.990 


0.943 


1 


1 


1 


1 


1 


1 


1 


1 


Tjy: Rejected i?^ 


1 


1 


1 


0.993 


0.330 


0.057 


0.007 


0.003 








ADF: Rejected H'^ 


1 


1 


1 


1 


1 


0.860 


0.687 


0.273 


0.043 


0.010 


PP : Rejected H^ 


1 


1 


1 


1 


1 


0.953 


0.857 


0.357 


0.027 





V/S : Accepted Ho 


1 


1 


0.943 


0.413 


0.080 


0.023 


0.007 












N = 5000 






















ARFIMA(l,d,l) 


d = -0.4 


d= -0.2 


d = 


ci = 0.2 


ci = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 


4> = -0.3 ; e = 0.7 






















Sm- Accepted Ho 


1 


1 


1 


1 


0.993 


0.703 


0.133 











KPSS: Accepted Hq 


1 


1 


0.937 


0.380 


0.037 


0.003 














LMC: Accepted Hq 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


T^: Rejected i?^ 


1 


1 


1 


1 


0.670 


0.043 














ADF: Rejected i?^ 


1 


1 


1 


1 


1 


0.990 


0.933 


0.460 


0.057 


0.010 


PP : Rejected //^ 


1 


1 


1 


1 


1 


1 


0.993 


0.607 


0.057 





V/S : Accepted Ho 


1 


1 


0.943 


0.200 





















Table 7: Comparisons of stationarity and nonstationarity tests from 300 independent replications of ARFIMA(1, d, 1) 

processes (with <f> = —0.3 and 9 = 0.7) for several values of d and A'^. The accuracy of tests is measured by the 
frequencies of trajectories "accepted as stationary" (accepted Hq or rejected Hq) among the 300 replications which 
should be close to 1 for d € (—0.5, 0.5) and close to for d e [0.5, 1.2] 



25 



iV = 200 


d = -0.4 


d= -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 
























S]\f- Accepted Ho 


1 


1 


1 


1 


1 


1 


0.977 


0.633 


0.093 


0.017 


KPSS: Accepted Hq 


1 


1 


0.997 


0.827 


0.460 


0.330 


0.220 


0.077 


0.017 


0.007 


LMC; Accepted Ho 


1 


1 


1 


0.837 


0.153 


0.103 


0.027 








0.317 


Tjv; Rejected H[^ 


1 


1 


0.997 


0.960 


0.720 


0.530 


0.173 


0.013 








ADF: Rejected 


1 


1 


1 


1 


1 


0.723 


0.620 


0.463 


0.280 


0.120 


PP : Rejected 


1 


1 


1 


1 


1 


0.997 


0.990 


0.880 


0.493 


0.113 


V/S : Accepted Ho 


1 


1 


0.99,3 


0.827 


0.427 


0.290 


0.187 


0.067 


0.027 


0.003 



N = 500 


d = -0.4 


d= -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 
























Spf: Accepted Hq 


1 


1 


1 


1 


1 


1 


0.947 


0.263 


0.017 





KPSS: Accepted Ho 


1 


1 


0.970 


0.693 


0.310 


0.130 


0.050 


0.013 


0.003 





LMC: Accepted Hq 


1 


1 


1 


0.670 


0.037 


0.003 











0.270 


Tn: Rejected i?^ 


1 


1 


1 


0.987 


0.803 


0.533 


0.207 


0.003 








ADF: Rejected H[^ 


1 


1 


1 


1 


1 


0.810 


0.820 


0.593 


0.300 


0.093 


PP : Rejected H'^ 


1 


1 


1 


1 


1 


1 


0.997 


0.913 


0.443 


0.057 


V/S : Accepted Hq 


1 


1 


0.990 


0.653 


0.207 


0.090 


0.053 


0.007 









N = 1000 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


d = 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d= 1.0 


d= 1.2 
























S]\r- Accepted Ho 


1 




1 


1 


1 


0.993 


0.917 


0.067 


0.003 





KPSS: Accepted Ho 


1 




0.990 


0.597 


0.157 


0.073 


0.023 


0.007 








LMC: Accepted Hq 


1 






0.553 

















0.177 


Tn: Rejected H^ 


1 






1 


0.910 


0.667 


0.213 











ADF: Rejected H'^ 


1 






1 


1 


0.900 


0.870 


0.630 


0.270 


0.087 


PP : Rejected H^^ 


1 






1 


1 


1 


1 


0.930 


0.427 


0.053 


V/S : Accepted Hq 


1 






0.580 


0.090 


0.027 


0.007 












iV = 5000 


d = -0.4 


d = -0.2 


d = 


d = 0.2 


(i= 0.4 


d = 0.5 


d = 0.6 


d = 0.8 


d = 1.0 


d = 1.2 


Sff: Accepted Ho 
KPSS: Accepted Ho 
LMC: Accepted Ho 


1 
1 
1 


1 
1 
1 


1 

0.977 
1 


1 

0.417 
0.213 


1 

0.050 



0.993 




0.653 



















Tjv: Rejected H[^ 
ADF: Rejected H^ 
PP : Rejected H^ 


1 
1 
1 


1 
1 
1 


1 
1 
1 


1 
1 
1 


0.977 
1 
1 


0.730 
1 
1 


0.073 
0.953 
1 




0.800 
0.983 




0.203 
0.320 




0.027 
0.007 


V/S : Accepted Ho 


1 


1 


0.987 


0.253 


0.007 


















Table 8: Comparisons of stationarity and nonstationarity tests from 300 independent replications of X'^'*'^-' processes 
for several values of d and A'^. The accuracy of tests is measured by the frequencies of trajectories "accepted as 
stationary" (accepted Ho or rejected Hq) among the 300 replications which should be close to 1 for d G (—0.5,0.5) 
and close to for d e [0.5, 1.2] 
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r=(USDl vs Deutsh-Mark Exchange Rate Return) 



Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 




dMS 


dROB 


^WAV 


[1 : 5963] 


S 


S 


S 


S 


s 


S 


SM 


5.5 


-0.2 


-0.031 


0.059 


0.057 


-0.007 


[5965 : 6313] 


S 


S 


s 


s 


s 


s 


SM 


3.4 


0.1 


0.034 


0.169 


0.122 


-0.015 


[6315 : 7173] 


S 


S 


s 


s 


s 


s 


SM 


5.3 


-0.4 


0.098 


0.140 


0.043 


0.019 


|r| = afts(USDl vs Deutsh-Mark Exchange Rate Return) 


Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 


d'R 


dMS 


dROB 


^WAV 


[1 ; 5963] 


S 


S 


S 


s 


NS 


NS 


LM 


9.5 


1.8 


0.294 


0.301 


0.344 


0.275 


[5965 : 6313] 


S 


S 


NS 


s 


NS 


NS 


LM 


3.6 


1.1 


-0.121 


0.153 


0.414 


-0.038 


[6315 : 7173] 


S 


S 


S 


s 


NS 


NS 


LM 


9.2 


1.8 


0.168 


0.417 


0.389 


0.410 


= (USDl vs Deutsh-Mark Exchange Rate Return)-^ 


Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 


d'R 


dMS 


dROB 


^WAV 


[1 : 5963] 


S 


S 


S 


s 


S 


S 


LM 


289.5 


10.7 


0.081 


0.258 


0.298 


0.078 


[5965 : 6313] 


S 


S 


S 


s 


NS 


S 


LM 


8.7 


2.3 


-0.018 


0.127 


0.431 


-0.096 


[6315 : 7173] 


S 


s 


S 


s 


NS 


S 


LM 


81.3 


7.1 


0.035 


0.411 


0.336 


0.428 



|r|* = (afes(USDl vs Deutsh-Mark Exchange Rate Return)) 





(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


ArgMaxe(S(|ri|»)) 


Sn 


Tn 


ADF 


PP 


KPSS 


LMC 


V/S 


K 


J 


d'R 




^ROB 


aWAv 




S 


S 


NS 


s 


NS 


NS 


SM 


3.5 


-0.5 


0.321* 


0.251 


0.256 


0.343 


ei^''>= 0.97 


S 


S 


S 


s 


NS 


NS 


SM 


8.7 


1.1 


0.293 


0.301* 


0.343 


0.275 




S 


s 


s 


s 


NS 


NS 


SM 


13.7 


2.3 


0.302 


0.300 


0.345* 


0.273 


^WAV)^^^^^ 


S 


s 


s 


s 


NS 


NS 


SM 


5.1 


1.1 


0.273 


0.298 


0.335 


0.379* 


0^2"*' =0.05 


S 


s 


NS 


NS 


S 


NS 


LM 


27.9 


-5.0 


0.246* 


0.078 


-0.005 


0.103 


e<"^)= 1.31 


S 


s 


S 


s 


NS 


NS 


LM 


4.9 


1.5 


-0.103 


0.166* 


0.446 


-0.072 


e</°^'=i.50 


S 


s 


S 


s 


NS 


S 


LM 


5.8 


1.8 


-0.092 


0.162 


0.450* 


-0.082 


§<'^^^>=0.03 


S 


s 


NS 


NS 


S 


NS 


LM 


30.9 


-5.4 


0.239 


0.113 


-0.030 


0.211* 


9</^)=0.63 


S 


s 


NS 


NS 


NS 


NS 


LM 


3.8 


0.7 


0.244* 


0.354 


0.333 


0.097 


§^*^«)= 1.44 


S 


s 


S 


S 


NS 


NS 


LM 


27.1 


3.7 


0.159 


0.436* 


0.387 


0.441 


e<'«°^'= 1.19 


S 


s 


S 


S 


NS 


NS 


LM 


14.9 


2.6 


0.168 


0.430 


0.394* 


0.430 


el^^v)^ 2.90 


S 


s 


S 


S 


NS 


S 


LM 


223.4 


13.3 


0.053 


0.291 


0.233 


0.475* 



Tabic 9: Results of stationarity, nonstationarity and V/S tests and the 4 long memory parameter estimators applied 
to several functionals / of USDl vs Deutsh-Mark Exchange Rate Return: from the top to bottom, f{x) = x, 
f{x) = \x\, f{x) = x^ and f{x) = with 6 maximizing the 4 different long memory parameter estimators ("S"' for 
"stationarity" decision and "NS" for "nonstationarity" decision). Statistics are applied to the 3 estimated stages of 
each trajectory (obtained from a change detection algorithm). 
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r=Dow Jones Transportation Index Return 



Segments 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skcwncss 


d 


Breaks 




Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 




dMS 


dROB 


flWAV 


[1 : 1271] 


S 


S 


S 


s 


s 


NS 


SM 


4.6 


0.0 


0.218 


0.174 


0.098 


0.198 


[1273 : 8531] 


S 


S 


s 


s 


s 


S 


SM 


21.7 


-0.8 


0.053 


0.002 


0.008 


-0.404 


[8533 : 12071] 


s 


S 


s 


s 


s 


S 


SM 


8.3 


-0.3 


0.002 


-0.015 


-0.034 


-0.038 


|r| = abs(Dow Jones Transportation Index Return) 


Segments 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skcwncss 


d 


Brealcs 




Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 


diR 


dMS 


flROB 


llWAV 


[1 : 1272] 


S 


S 


S 


s 


NS 


NS 


LM 


6.1 


1.5 


0.154 


0.320 


0.270 


0.166 


[1273 : 8532] 


S 


S 


s 


s 


NS 


NS 


LM 


57.3 


4.4 


0.322 


0.260 


0.240 


0.168 


[8533 : 12071] 


S 


S 


s 


s 


NS 


NS 


LM 


16.3 


2.5 


0.405 


0.476 


0.496 


0.374 


r-^ = (Dow Jones Transportation Index Return)^ 


Segments 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skcwncss 


d 


Breaks 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 


d'R 


dMS 


dROB 


fjWAV 


[1 : 1272] 


S 


S 


S 


s 


NS 


S 


LM 


32.3 


4.4 


0.158 


0.284 


0.231 


0.231 


[1273 : 8532] 


S 


S 


s 


s 


S 


NS 


LM 


2301.5 


39.9 


0.334 


0.122 


0.093 


0.118 


[8533 : 12071] 


S 


NS 


s 


s 


NS 


NS 


LM 


459.0 


15.5 


0.416 


0.452 


0.434 


0.356 



|r|* = (a6s(Dow Jones Transportation Index Return))* 





(Non)Stationarity Test 


LRD 


Kurtosis 


Skcwncss 


d 


ArgMaX(,(d(|ri|*) 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 


d'R 


dMS 


dROB 


dWAV 


e^/^>= 1.83 


S 


S 


S 


s 


NS 


S 


LM 


25.0 


3.8 


0.252* 


0.291 


0.237 


0.202 


0.45 


S 


S 


NS 


s 


NS 


NS 


LM 


2.7 


0.4 


0.118 


0.331* 


0.290 


0.047 


q(ROB)^ 0.36 

e^^^■A^)= 0.03 


S 
S 


S 
S 


NS 
NS 


s 
s 


NS 
NS 


NS 
NS 


LM 
LM 


2.9 
12.8 


-0.4 
-3.4 


0.118 
0.149 


0.329 
0.257 


0.291* 

0.260 


0.237 
0.327* 


6^^-"'= 2.06 


S 


S 


S 


s 


S 


NS 


LM 


2551.6 


42.6 


0.355* 


0.113 


0.086 


0.110 


e</^^)= 0.68 


S 


S 


S 


s 


NS 


NS 


LM 


10.6 


1.6 


0.308 


0.276* 


0.261 


0.135 


gtflOB)^ 0.65 
9^^'^^^= 1.29 


S 
S 


S 
S 


S 
S 


s 
s 


NS 
NS 


NS 
NS 


LM 
LM 


9.2 

246.8 


1.4 
10.1 


0.303 
0.330 


0.276 
0.227 


0.261* 

0.200 


0.129 
0.504* 


e^^^'= 0.66 


S 


NS 


S 


s 


NS 


NS 


LM 


5.4 


1.1 


0.444* 


0.435 


0.461 


0.374 


d\t''^= 1.38 


s 


S 


S 


s 


NS 


NS 


LM 


64.8 


5.2 


0.402 


0.492* 


0.499 


0.391 


q{ROB)^ 1.22 


s 


S 


S 


s 


NS 


NS 


LM 


36.2 


3.8 


0.400 


0.489 


0.502* 


0.387 


e<3^-"^)= 2.75 


s 


s 


s 


s 


NS 


NS 


LM 


1698.7 


35.8 


0.407 


0.315 


0.287 


0.466* 



Tabic 10: Results of stationarity, nonstationarity and V/S tests and the 4 long memory parameter estimators 
applied to several functionals / of DowJones Transportation Index Return: from the top to bottom, f(x) = x, 
f{x) = |a;|, f{x) = x^ and f{x) = \xf with 6 maximizing the 4 different long memory parameter estimators ("S"' for 
"stationarity" decision and "NS" for "nonstationarity" decision). Statistics are applied to the 3 estimated stages of 
each trajectory (obtained from a change detection algorithm). 
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r=Dow Jones Utilities Index Return 



Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sjv 


Tjv 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 






^ROB 


dWAV 


[1 : 1152] 


S 


S 


S 


s 


S 


S 


SM 


7.3 


0.6 


0.191 


0.037 


-0.132 


0.222 


[1153 : 8748] 


S 


s 


S 


s 


S 


s 


SM 


43.2 


-1.3 


0.094 


0.025 


0.001 


0.043 


[8749 : 12071] 


S 


s 


s 


s 


s 


s 


SM 


13.0 


0.0 


0.026 


0.024 


0.001 


-0.032 


|r| = abs(Dow Jones Utilities Index Return) 


Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 




Tjv 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 




dMS 


d^OB 


flWAV 


[1 ; 1152] 


S 


S 


s 


s 


NS 


NS 


LM 


11.9 


2.4 


0.283 


0.287 


0.316 


0.225 


[1153 : 8748] 


s 


s 


s 


s 


NS 


NS 


LM 


127.4 


5.9 


0.134 


0.301 


0.304 


0.184 


[8749 : 12071] 


s 


s 


s 


s 


NS 


NS 


LM 


25.5 


3.4 


0.417 


0.559 


0.484 


0.595 


= (Dow Jones Utilities Index Return)'^ 


Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sjv 


Tjv 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 




dMS 


d^OB 


fjWAV 


[1 : 1152] 


S 


S 


s 


s 


NS 


S 


LM 


63.1 


6.7 


0.250 


0.253 


0.212 


0.270 


[1153 : 8748] 


S 


s 


s 


s 


NS 


NS 


LM 


5322.4 


67.8 


0.130 


0.100 


0.100 


0.100 


[8719 : 12071j 


s 




s 


s 


NS 


NS 


LM 


281). (i 


ii.O 


o.rjio 


0. 1()8 


0. 123 


0.513 



[r]" = (afcs(Dow Jones Utilities Index Return)) 





(Non)Stationarit,y Tost 


LRD 


Kurtosis 


Skewness 


d 


ArgMa.x,(,7(,r,,")) 




T.x 


ADF 


I'l' 


Kl'SS 


LMC 


I '/S 












,IU-AV 


e['"-^ = 0.39 


s 


S 


NS 


s 


NS 


NS 


LM 


3.3 


-0.1 


0.354* 


0.262 


0.327 


0.145 


e[^''> = 1.09 


s 


S 


S 


s 


NS 


NS 


LM 


14.6 


8.4 


0.215 


0.288* 


0.308 


0.234 




s 


S 


NS 


s 


NS 


NS 


LM 


4.5 


0.8 


0.311 


0.276 


0.336* 


0.396 


§<"'^^) = 0.63 


s 


s 


S 


s 


NS 


NS 


LM 


4.8 


1.0 


0.310 


0.278 


0.336 


0.398* 


= 3.00 


s 


s 


s 


s 


s 


S 


LM 


7320.6 


84.9 


0.165* 


0.015 


0.017 


0.040 


= 0.61 


s 


s 


s 


s 


NS 


NS 


LM 


9.1 


1.2 


0.113 


0.330* 


0.327 


0.113 




s 


s 


s 


s 


NS 


NS 


LM 


13.0 


1.6 


0.117 


0.330 


0.327* 


0.113 


e<^^^) = 1.84 


s 


s 


s 


s 


NS 


NS 


LM 


4386.6 


59.1 


0.125 


0.130 


0.129 


0.377* 


9^^^* = 2.69 


s 


NS 


s 


s 


NS 


NS 


LM 


683.1 


22.5 


0.527* 


0.394 


0.344 


0.426 


= 0.95 


s 


s 


s 


s 


NS 


NS 


LM 


21.6 


3.1 


0.415 


0.560* 


0.483 


0.544 


= 1.10 


s 


s 


s 


s 


NS 


NS 


LM 


35.2 


4.2 


0.421 


0.557 


0.485* 


0.364 


e<3'^'*^) = 1.03 


s 


s 


s 


s 


NS 


NS 


LM 


28.2 


3.7 


0.419 


0.559 


0.484 


0.723* 



Tabic 11: Results of stationarity, nonstationarity and V/S tests and the 4 long memory parameter estimators applied 
to several functionals / of Dow Jones Utilities Index Return; from the top to bottom, f{x) = x, f{x) = \x\, f{x) = x^ 
and f{x) = \xf with maximizing the 4 different long memory parameter estimators ("S"' for "stationarity" decision 
and "NS" for "nonstationarity" decision). Statistics are applied to the 3 estimated stages of each trajectory (obtained 
from a change detection algorithm). 
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r=Nasdaq Industrials Index Return 



Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skcwncss 


d 




Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 




dMS 


dROB 


fjWAV 


[1 : 7160] 


S 


S 


S 


s 


s 


S 


SM 


20.7 


-1.5 


0.141 


0.073 


0.092 


-0.202 


[7161 : 8320] 


S 


S 


S 


s 


s 


S 


SM 


4.6 


0.0 


0.012 


0.070 


0.116 


0.014 


[8321 : 10480] 


s 


S 


s 


s 


NS 


S 


SM 


10.4 


-0.3 


0.045 


0.078 


0.082 


-0.045 


|r| = abs (Nasdaq Industrials Index Return) 


Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skcwncss 


d 




Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 


diR 


dMS 


flROB 


llWAV 


[1 : 7160] 


S 


S 


S 


s 


s 


NS 


SM 


52.4 


4.4 


0.361 


0.309 


0.287 


0.274 


[7161 : 8320] 


S 


S 


s 


s 


NS 


NS 


LM 


7.4 


1.6 


0.284 


0.532 


0.504 


0.385 


[8321 : 10480] 


S 


NS 


s 


NS 


NS 


NS 


LM 


18.3 


3.0 


0.516 


0.761 


0.606 


0.668 


= (Nasdaq Industrials Index Return)^ 


Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skcwncss 


d 


Sn 


Tn 


ADF 


PP 


KPSS 


LMC 


V/S 


K 


s 


d'R 


dMS 


dROB 


fjWAV 


[1 : 7160] 


S 


S 


S 


s 


s 


NS 


SM 


1356.8 


31.4 


0.381 


0.146 


0.114 


0.100 


[7161 : 8320] 


S 


S 


s 


s 


NS 


NS 


LM 


49.0 


5.4 


0.304 


0.466 


0.378 


0.432 


[8321 : 10480] 


s 


NS 


s 


s 


NS 


NS 


LM 


140.0 


10.0 


0.498 


0.786 


0.544 


0.708 



|r|* = (a6s(Nasdaq Industrials Index Return)) 





(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


ArgMaX(,(5(|ri|»)) 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


5' 


d'R 




dROB 


dWAV 


e['"'> = 1.04 


S 


S 


S 


s 


s 


NS 


LM 


63.7 


4.9 


0.396* 


0.304 


0.281 


0.325 


ej*^*') = 0.67 


S 


S 


S 


s 


NS 


NS 


LM 


10.2 


0.7 


0.188 


0.329* 


0.325 


0.293 


^ROB) ^ g 


S 


S 


s 


s 


NS 


NS 


LM 


6.3 


0.9 


0.178 


0.326 


0.328* 


0.275 


e<^'^^) = 0.83 


S 


S 


s 


s 


NS 


NS 


LM 


22.4 


2.6 


0.199 


0.324 


0.311 


0.587* 


e^'-"' = 2.83 


s 


S 


s 


s 


NS 


S 


LM 


142.0 


9.8 


0.317* 


0.374 


0.276 


0.263 


ei"^^ = 1.03 


s 


S 


s 


s 


NS 


NS 


LM 


7.9 


1.7 


0.284 


0.532* 


0.501 


0.388 


q(ROB) ^ 


s 


S 


NS 


s 


NS 


NS 


LM 


4.2 


0.9 


0.300 


0.517 


0.517* 


0.340 


e^^^^^ = 1.87 


s 


S 


S 


s 


NS 


NS 


LM 


39.8 


4.8 


0.299 


0.479 


0.395 


0.432* 


e^'"-'' = 2.60 


s 


NS 


S 


s 


NS 


S 


LM 


256.9 


14.3 


0.548* 


0.272 


0.479 


0.669 


= 1.70 


s 


NS 


S 


s 


NS 


NS 


LM 


89.1 


7.7 


0.504 


0.801* 


0.575 


0.739 


= 1.13 


s 


NS 


s 


s 


NS 


NS 


LM 


26.0 


3.7 


0.526 


0.772 


0.608* 


0.671 


4^^^^) = 1.26 


s 


NS 


s 


s 


NS 


NS 


LM 


36.0 


4.5 


0.532 


0.782 


0.606 


0.760* 



Table 12: Results of stationarity, nonstationarity and V/S tests and the 4 long memory parameter estimators applied 
to several functionals / of Nasdaq Industrials Index Return: from the top to bottom, f{x) = x, f{x) = \x\, f{x) = x^ 
and f{x) = I a; I* with maximizing the 4 different long memory parameter estimators ("S"' for "stationarity" decision 
and "NS" for "nonstationarity" decision). Statistics are applied to the 3 estimated stages of each trajectory (obtained 
from a change detection algorithm). 
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r=Nikkei 225A Index Return 



Segments 
Breaks 


(Non)Statioiiarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sn 


Tm 


ADF 


pp 


KPSS 


LMC 


v/s 


K 








dROB 


dWAV 


[1 : 6672] 


S 


S 


S 


s 


s 


s 


SM 


12.6 


-0.6 


0.083 


0.067 


0.084 


0.022 


[6673 : 10400] 


S 


S 


S 


s 


s 


s 


SM 


63.7 


-2.3 


-0.021 


-0.016 


-0.013 


-0.039 


[10401 : 15919] 


S 


S 


s 


s 


s 


s 


SM 


9.0 


-0.1 


0.033 


0.047 


-0.005 


-0.015 



\r\ = af)s(Nikkei 225A Index Return) 



Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 




s 


d'R 


dMS 


^ROB 


^WAV 


[1 : 6672] 


S 


S 


S 


s 


NS 


NS 


LM 


26.1 


3.3 


0.302 


0.343 


0.313 


0.218 


[6673 : 10400] 


S 


S 


s 


s 


NS 


NS 


LM 


150.9 


7.5 


0.196 


0.346 


0.304 


0.321 


[10401 : 15919] 


S 


s 


s 


s 


NS 


NS 


LM 


17.0 


2.6 


0.413 


0.415 


0.431 


0.335 



r2 = (Nikkei 225A Index Return)^ 



Segments 
Breaks 


(Non)Stationarity Test 


LRD 


Kurtosis 


Skewness 


d 


Sn 


Tn 


ADF 


pp 


KPSS 


LMC 


V/S 


K 


s 




dMS 


d^OB 


^WAV 


[1 : 6672] 


S 


S 


S 


s 


NS 


NS 


LM 


427.8 


16.9 


0.275 


0.241 


0.267 


0.080 


[6673 : 10100] 


S 


s 


s 


s 


S 


NS 


LM 


2610.8 


48.0 


0.2.30 


0.116 


0.1.54 


0.117 


[10101 : 15919j 


s 


s 


s 


s 


NS 


NS 


LM 


2:-!5.."') 


12.:-! 


0.:-!81 


0.:-!IJ() 


0. :-!():-! 


0.377 



jrl" = (afe6'(Nikkei 225A Index Return))' 





(Non)Stat,ionarit,y Tost 


LRD 


Kurtosis 


Skewness 


d 


ArgMa.x„(,7(,r,,")) 




T.x 


ADF 


I'l' 


Kl'SS 


LMC 


I '/S 














e*'-"' = 1.53 


s 


S 


S 


S 


NS 


NS 


LM 


149.8 


9.0 


0.323* 


0.296 


0.296 


0.201 


= 0.86 


s 


S 


S 


S 


NS 


NS 


LM 


15.4 


4.4 


0.286 


0.345* 


0.311 


0.213 


^ROB) ^ ^ 


s 


S 


S 


S 


NS 


NS 


LM 


27.1 


3.4 


0.303 


0.342 


0.313* 


0.218 


e^^^^^ = 1.30 


s 


s 


S 


S 


NS 


NS 


LM 


74.8 


6.1 


0.273 


0.322 


0.307 


0.622* 


= 3.00 


s 


s 


S 


S 


S 


S 


LM 


3487.4 


58.3 


0.252* 


0.037 


0.042 


0.045 


Pf^^ = 0.84 


s 


s 


S 


S 


NS 


NS 


SM 


55.5 


4.1 


0.180 


0.353* 


0.304 


0.154 


0(ROB) ^ Q 


s 


s 


S 


S 


NS 


NS 


SM 


87.1 


5.3 


0.186 


0.352 


0.305* 


0.035 


e<^^^) = 1.64 


s 


s 


S 


S 


S 


NS 


SM 


1697.1 


35.8 


0.221 


0.220 


0.222 


0.465* 


e\'^^ = 1.23 


s 


NS 


S 


S 


NS 


NS 


LM 


35.7 


4.1 


0.467* 


0.412 


0.426 


0.386 


= 0.87 


s 


s 


S 


S 


NS 


NS 


LM 


11.0 


2.0 


0.429 


0.415* 


0.428 


0.351 


giROB) ^ 


s 


s 


S 


S 


NS 


NS 


LM 


17.0 


2.6 


0.413 


0.415 


0.431* 


0.335 


9(WAV) ^ ^ 


s 


NS 


S 


S 


NS 


NS 


LM 


40.3 


4.4 


0.467 


0.411 


0.424 


0.425* 



Tabic 13: Results of stationarity, nonstationarity and V/S tests and the 4 long memory parameter estimators applied 
to several functional / of Nikkei 225A Index Return: from the top to bottom, f{x) = x, f{x) = \x\, f{x) = x^ and 
f{x) = \xf with 9 maximizing the 4 different long memory parameter estimators ("S"' for "stationarity" decision and 
"NS" for "nonstationarity" decision). Statistics are applied to the 3 estimated stages of each trajectory (obtained 
from a change detection algorithm). 
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